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ABSTRACT 


A  simplified  two-degree-of-f reedom  simulation  in 
heave  and  pitch  of  the  XR-3  Captured  Air  Bubble  (CAB) 
testcraft  is  developed  first  in  a  non- linear  form  and 
is  subsequently  linearized  about  the  steady-state 
operating  point.  The  model  is  validated  against  the 
six-degree-of-f reedom  XR-3  Loads  and  Motions  program 
in  the  time  domain;  the  linear  system  is  then 
transformed  into  the  frequency  domain  by  complex 
matrix  inversion.  A  comparison  is  made  between  the 
linear  system  frequency  response  predicted  by  the 
two-degree-of- freedom  model  and  that  of  the  Loads  and 
Motions  program,  as  well  as  actual  testcraft  data, 
reduced  to  the  frequency  domain  by  a  Fast  Fourier 
Transform       (FFT)       technique.  This        comparison        of 

frequency  response  curves  highlights  a  non-linear  mode 
of  CAB  craft  behavior  of  possible  interest  to 
habitability   and   seaworthiness   studies. 
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LIST  OF  ABBREVIATIONS 


PHYSICAL  CONSTANTS 


Gravitational  Constant 


(32. 2  ft/sec  ) 


2 .   Gamaa 


Satio  of  Specific  Heats  of  Air 

( 1 . 4, Dim en si onl ess) 


3 .   Rho 


4.   Rhoa 


Density  of  Water 


(2.0  Slugs/ft  ) 


Density  of  Air 


(.002378  Slugs/ft  ) 
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CRAFT    PARAMETERS 


1.   Ab 


2.   Al 


3.   As 


Plenum  Upper  Surface  Area 


(200  ft  ) 


Leakage  Area, Stern  Seal 


Sidewall  Area 


(.438  ft   ) 


(13.75  ft  ) 


4.   Cf 


Sidewall  Skin  Friction  Coefficient 

( .  95  ,Dioensionless) 


5.      Ill 


Y-Axis   Moment   of    Inertia 


(9620    lbm-ft    ) 


6.      L 


Craft  Waterline    Length 


(20.0  f-c) 


7.   N 


Number  of  Lift  Fans 


(5) 


8.   £T_c 


Static  Fan  Flow  Rate 


(35  ft/sec) 


9.   Ws 


Sidewall  Width 


(.9375  ft) 
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10.       Width  Craft   Beam    Width 

—  (10.0    ft) 


11.       Vn  Empty   Plenum    Volume 

—  3 

(333  ft  ) 


12.   Zs  Keel  to  Center  of  Gravity  Height 

(2.5  ft) 


13.   M  Mass  of  Craft 

(lb  in) 


14.   Mf C  Residual  Moment  about  Y-Axis 

(ft-lbf) 


15.      Xcg  Location      of    Center    of   Gravity    from 

Craft  Geometric    Center 

(ft) 


16.   Xcp  Location    of    Craft   Center   of 

Pressure    from    Craft.    Geometric 
Center 

(ft) 
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SEA-STATE    VARIABLES 


1.   C 


Absolute  Wave  Velocity 


(ft/sec) 


2.   Lambda 


Incident  Wavelength 


(ft) 


3.  v 


4.   v 


Wave  Encounter  Frequency 


Incident  Wave  Frequency 


(rad/sec) 


(rad/sec) 


5.   Wvht 


Wave  Height 


(ft) 
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STATE  VARIABLES 


1 .   Mb  Plenum  Air  Mass 

(lbm) 


2.   Alfa  Angular  Acceleration  of  Craft  about 

Y-Axis 


2 
(rad/sec  ) 


3.   Mbdot  Mass  Flow  Rate  into/cut  of  Plenum 

(lbm/sec) 


4.   Omega  Radial    Velocity   of   Craft   about 

~  Y-Axis 

(rad/sec) 


5.   Pbbar  Plenum  Gauge  Pressure 

2 
(lbf/ft  ) 


6.   Theta  Pitch  Angle  about  Y-Axis 

(rad/sec) 


7.   Zdot  Vertical    Velocity   of   Center   of 

Gravity 

(ft/sec) 


8.   Z2dot  Vertical   Acceleration  of  Center  of 

Gravity 

2 
(ft/sec  ) 


9.   Z  Vertical   Displacement  of  Center  of 

Gravity  above  Water  Line 

(ft) 
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SYSTEM  VARIABLES 


1.   Aseall 


Wetted  Area  of  Bowseal 


2 
(ft  ) 


2.   Aseal2 


Wetted  Area  of  Sternseal 


(ft) 


3 .  Hbuoy 

4 .  Hbuo^j 


Total  Buoyant  Force 


(lbf) 


Buoyant  Force  due  to  immersed   Area 
forward  of  Center  of  Gravity 

(lbf) 


5 .   Hbuoy2 


Buoyant  Force  due  to  Immersed   Area 
aft" of  Center  of  Gravity 

(lbf) 


6 .   Hbuoyw 


Additional   Force    due    to    Wave 
Buoyancy 

(lbf) 


7 .  Hdamp 

8.  Hdaapl 


Total  Viscous   Damping   Force 


(lbf) 


Damping  Force  due  to  Viscous 
Damping  forward  of  Center  of 
Gravity 

(lbf) 


9 .   Hdaap2 


Damping   Force    due   to    Viscous 
Damping  aft  of  Center  of  Gravity 

(lbf) 
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10.   Hplan 


Vertical  Planinq   Force    due   to 


vertical  Planing 
Hydrodynamic  Lift 


(lbf) 


1 1 .   Hpres 


Vertical  Pressure  Force,  acting   at 
Center  of  Pressure 

(lbf) 


12.   Hseall 


Vertical  Force  due  to  Bowseal 


(lbf) 


13.   Hseal2 


Vertical  Force  due  to  Stsrnseal 


(lbf) 


14.       Ld 


Draft      , measured        at        Center        of 
Gravity 

(ft) 


15.   Ld1 


Draft,  measured  at  bow 


(ft) 


16.   Ld2 


Draft,  measured  at  stern 


(ft) 


17.   Ldbarl 


Average  Draft  forward  of  Center   of 
Gravity 

(ft) 


18.   Ldbar2 


Average   Draft   aft   of   Center   of 
Gravity 

(ft) 


1  9 .      PbpaC 


Initial   Plenum   Pressure  Ratio 

(dime  nsionl  ess) 


20.      Pbuoy 


Total   Buoyant   Pitch    rtoment 


(f t-lbf) 


2 1 .   Pbuoyw 


Buoyant  Pitch  Moment  due  to  Haves 

(f  t-lbf) 
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22.   Pbuoyl 


Buoyant  Pitch  Moment  due  to 
Immersed  Area  forward  of  Center  of 
Gravity 

(ft-lbf) 


23 .   Pbuoy2 


Buoyant  Pitch  Moment  due  to 
Immersed  Area  aft  of  Center  of 
Gravity 

(ft-lbf) 


2<f. 


25. 


PdamjD 
Pdampl 


Total  Viscous  Pitch  Damping  Moment 

(ft-lbf) 


Viscous  Pitch  Damping  Moment  due  to 
Area  forward  of  Center  of  Gravity 

(f t-lDf) 


26.   Pdamg2 


Viscous  Pitch  Damping  Moment  due  ~o 
Area  aft  of  Center  of  Gravity 

(ft-lbf) 


27.   P£lan 


Planing  Moment  due  to   Hydrcdynamic 
Lift 

(ft-lbf) 


28.   Ppres 


Pitch  Pressure  Moment 


(ft-lbf) 


29.   Pseall 


Pitch  moment  due  to  Bowseal 


(f t-lof) 


30.   Pseal2 


Pitch  moment  due  to  Sternseal 

(ft-lbf) 


3  1 .   Qout 


32.   Qin 


Volume  Flow  out  of  Plenum 


(ft/sec) 


Volume  Flow  into  the  Plenum 


(ft/sec) 
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33.   Vb  Net  Plenum  Volume 

3 
(ft  ) 


34.   Vbar  1  Average  Vertical   Velocity   forward 

of  Center  of  Gravity 

(ft/sec) 


35.   Vbar2  Average  Vertical   Velocity   aft   of 

Center   of   Gravity 

(ft/sec) 


36.   Xseal  Length  of  Bowseal  in   Contact   with 

the  Water  Surface 

(ft) 
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INFLUENCE  COEFFICIENTS 


1.   Dhbth  Partial   derivative   of    vertical 

buoyant  force  with  respect  to  Theta 

(  lbf) 


2.   Dhbz  Partial   derivative   of    vertical 

buoyant  force  with  respect  to  z 

(lbf/ft) 


3.   Dhglth  Partial   derivative   of    vertical 

force   due   to  planing  with  respect 
to  Theta 

(lbf) 


4.   Dhjgth  Partial    derivative    of   vertical 

pressure   force   with   respect   to 
Theta 

(  lbf) 


5.   Dhpmb  Partial   derivative   of    vertical 

pressure  force  with  respect  tc  Mb 

(    lbf/lbm) 


6.      DhjDZ  Partial      derivative        of        vertical 

pressure    force    with    respect    to    z 

(lbf/ft) 


7.      Dhsz 1  Partial        derivative        of        flcwseal 

vertical    force    with   respect    tc   z 

(  lbf/ft) 


8.   Dhsth.1  Partial    derivative   of    Bowseal 

vertical    force   with   respect   to 
Theta 

(lbf) 
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9.   Dhsab  Partial   derivative   of    Bowseal 

vertical  force  with  respect  to  Mb 

(lbf/lbm) 


10.   Dhsz2  Partial    derivative   of   Sternseal 

vertical  force  with  respect  tc  z 

(  ltf/ft) 


11.   ffhsth2  Partial   derivative   of   Sternseal 

vertical   force   with   raspect   to 
Theta 

(lbf) 


12.   Dhsz  Sum  of  Dhszl  and  Dhsz2 

(lbf/ft) 


13.       Dhsth  Sum    of   Dhsthl    and    Dhsth2 

(lbf) 


14.   Dmbmb  Partial  derivative   of   Plenum   Air 

Mass  Flow  Rate  with  respect  tc  Mb 

-1 
(sec) 


15.   Dmbth  Partial  derivative   of   Plenum   Air 

Mass   Flow   Rate   with   respect   to 
Theta 

(  lbm/sec) 


16.   Dmbz  Partial   derivative   of   Plenum  Air 

Mass  Flow  Rate  with  respect  tc  z 

(lbm/f  t-sec) 


17.      Dpbab  Partial      derivative      of      Pbbar   with 

respect    to   Plenum    Air   Mass,    Mfc 

2 
(  lbf/lbm-ft) 


18.   Dpj2th  Partial   derivative   of   Pbbar  with 

respect  to  Theta 

2 
(ltf/ft) 
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19.   Dpbz  Partial   derivative   of   Pbbar  with 

respect  to  z 

3 
(lbf/ft  ) 


20.   Dobjth  Partial  derivative  of  buoyant  pitch 

moment  with  respect  to  Theta 

(ft-lbf) 


21*   Dpbyz  Partial  derivative  cf  buoyant  pitch 

moment  with  respect  to  z 

(lbf) 


22.   Dpplth  Partial  derivative  cf  Planing  pitch 

moment  with  respect  to  Theta 

(ft-lbf) 


23.   Dppmb  Partial    derivative   of   pressure 

pitch  moment  with  respect  to  Mb 

(f  t-lbf/lbm) 


24.      DpjDth  Partial        derivative        of      pressure 

pitch    moment    with    respect    to    Theta 

(ft-lbf) 


25.      Dppz  Partial        derivative        of      pressure 

pitch    moment    with    respect    to    z 

(lbf) 


26.   Dpsab  Partial  derivative  of  Bowseal  pitch 

moment  with  respect  to  Ab 

(  ft-lbf/lbin) 


27.   DpsthJ  Partial  derivative  of  Bowseal  Ditch 

moment  with  resoect  to  Tneta 

(  ft-lbf) 


28.      Dgsth2  Partial        derivative      of      Sternseal 

pitch    moment   with    respect    to    Theta 

(ft-lbf) 


29.   Dpsz  1  Partial  derivative  of  Bowseal  pitch 

moment  with  respect  to  z 

(  lbf) 
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30.   Ppsz2  Partial    derivative   of   Sternseal 

pitch  moment  with  respect  to  z 

(lbf) 


31.   Ppjnb  Partial   .  derivative         p 


acceleration  with  respect  to 

2 
(  1/(lbm-sec)) 


32.      Dpth  Partial    derivative    of    pitch 

acceleration  with  respect  to  Ineta 

-2 
(  sec) 


33.      Dgz  Partial  derivative  of        pitch 

acceleration    with    respect    to    z 

2 
(    1/(ft-sec)) 


34.   Dzmb  Partial   derivative   of   vertical 

acceleration  with  respect  to  Mb 

2 
(  ft/lbm-sec) 


35.   Dzth  Partial   derivative   of  of  vertical 

acceleration  with  resDect  to  Iheta 

2 
(ft/sec) 


36.   Dzz  Partial   derivative   of   vertical 

acceleration  with  respect  to  z 

-2 
(  sec) 
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I.   INTRODUCTION 


The  frequency  domain  provides  a  asaful  method  of  study 
for  systems  such  as  the  Captured  kir  Buoble  (CAB)  type 
surface  effect  ship,  wnich  are  excited  by  forces  *no^: 
statistics  of  amplitude  and  rrequency  may  be  well-xnown,  but 
wnose  time  domain  phase  relationships  ire  quite 
indeterminan t .  In  addition,  the  frequency  response  of  such 
hign-speed  ships  should  be  closely  examinee  for  resonances 
near  wave  encounter  frequencies  most  likely  to  be 
experienced.  Due  to  its  relatively  shallow  drift,  the  CA3 
craft  can  easiay  be  excited  into  a  noa-iiuear  regicn  of 
operation  in  which  sidewa.ll  gapping  or  wave/plenum  up]  - 
surface  contact  may  occur;  tnese  non-linearities  may  excite 
modes  of  operation  which  are  of  importance  to  habitability, 
and  may  limit  the  craft's  speed  during  sea-state   operation. 

Boota   in   ref£l7  analyzed  the  fr j^aency  domain  response 
of  the  X.i-3  CAt3  craft  utilizing  the  i,oads  and  Motion  program 
developed   by  Oceanics, Inc. ,   for   the   bell   1003  C!\3    ind 
modified  extensively  by  Leo  and  Boncal  (ref  £-7)  and   Kenzel 
(rer£jj)  ,   ror   the   X?-3.    However,   the   Loads  and  Koti  );: 
program  is  an  extremely   inefficient   tool   to   utilize   for 
Fourier   series   analysis;   not   only   ar i    waves  numerically 
integrated  along  the  hull  utilizing  a   table   look-up,   once 
for   eacn   frequency  component,  and  four  times  for  aach  time 
step,  out   the  higher  frequency   components   also   drive   the 
variable-ste  t-   integrators   to   extremely   small  step  siz 
The   net   effect   of   this   technique   is    the    exortit 
utilization    of   CPU   time   in   a   computer    .  ilaticn 
sea-stat3  conditions.   In  addition.,  th    time   constants   ot 
the   crait   lesponse   are   relatively   long       some  to  i  ., 


particularly  pitch;  therefore  still  further  CPU  time  is 
wasted  attaining  steady-state  conditions.  This  inefficiency 
of  CPU  time  has  been  noted  by  Mitchell  in  ref  Fu J ,  who 
documented  over  40  percent  of  CPU  time  utilized  by  the  XR-3 
Loads  and  Motion  program  to  the  WAVES  subroutine. 
Accordingly,  one  is  reduced  to  performing  a  point-by-point 
freguency  analysis  w;Lth  a  single  wave  freguency  at  a  time. 
Since,  under  these  conditions,  the  system  is  approximately 
linear,  and  superposition  is  valid  for  sufficiently  small 
wave  amplitudes,  the  utilization  of  a  linearized,  albeit 
highly  simplified,  program  presents  itself  as  a  most 
desirable  alternative.  The  effect  of  non-linearities  nay  be 
added  to  the  linear  model  with  a  describing  function,  and 
the  system  investigated  in  the  complex  plane  utilizing  the 
modified  Nyguist  criteria.  Finally,  a  simplified  model 
provides  a  starting  point  for  the  investigation  of 
sensitivity  of  craft  operation  to  the  variation  of  assorted 
parameters  through  the  technique  of  root  locus  plotting,  and 
can  be  the  heart  of  the  development  of  any  automatic  control 
system,  optimal  or  classical. 

Gerba  and  Thaler,  in  ref  \5j  ,  initiated  the  development 
of  the  linearized  XR-3  model  in  heave  only.  In  this  work, 
that  model  is  further  developed  to  include  the  pitch  mode, 
first  in  a  non-linear  form,  from  which  a  linearized  version 
is  derived.  The  wave  forcing  functions  acting  along  the 
hull  of  the  XR-3  are  developed  in  the  form  of  a  quasi-linear 
B  vector  in  state-space,  and  the  linear  system  is  then 
transformed  into  the  frequency  domain  by  complex  matrix 
inversion,  with  significant  savings  in  CPU  time.  This 
linear  frequency  response  is  then  compared  with  ref  fij  to 
provide  additional  insights  into  the  XR-3  system  behavior 
not  readily  apparent  in  the  more  complex  Loads  and  Motions 
program. 

A   technique   utilizing   a   mini-computer   and   spectral 
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analyzer  is  developed  for  the  fast  and  efficient  reduction 
of  recorded  XR-3  data  from  the  time  domain  to  the  freguency 
domain  utilizing  a  Fast  Fourier  Transform  (FFT) .  Actual 
craft  data  is  compared  with  the  predicted  linear  system 
frequency  response. 
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II.   NON^LINEAR  XR-3  HEAVE-PITCH  EQUATIONS  OF  MOTION 


This  chapter  develops  the  non-linear  equations  of 
motions  of  the  XR-3  in  heave  and  pitch,  following  the 
approach  utilized  by.  Gerba  and  Thaler  in  the  development  of 
the  simplified  heave- only  model  in  ref  £51  . 


A.   ASSUMPTIONS 

As  in  ref  [51  ,  the  following  assumptions  are  made: 

1.  constant  leakage  area; 

2.  equilibrium  speed  conditions; 

3.  plenum  pressure  and  volume  related  through  an  adiabatic 
process ; 

The   following   assumptions   and   limiting   conditions   were 
removed: 


4.  vertically   coincident    geometric 
pressure  and  center  of  gravity; 

5.  negligible  pitch  changes; 

6.  negligible  wave  action; 

7.  negligible  added  mass  effects; 


center , center 


and  the  following  assumptions  are  added: 
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8.  negligible  wave  action  within  plenum; 

9.  wave    action    due   to   head-on   or   following   seas 
(negligible  roll  or  yaw); 

In  addition  some  simplifying  assumptions  regarding  the  craft 
geometry  were  made.  Fig  1  shows  the  simplified  geometry 
of  the  craft  assumed  in  this  model.  Note  that  Center  of 
Gravity  (C.G. )  and  Center  of  Pressure  (C.P.)  are  arbitrarily 
located  and  not  to  scale. 


B.   COORDINATE  SYSTEM 


Pig  1  shows  the  coordinate  system  used  in  this  model. 
Specifically,  the  coordinate  system  is  attacaed  to  the 
geometric  center  of  the  craft  at  the  calm  watarline,  and 
oriented  in  a  right-hand  system  such  that  x  is  positive 
forward,  y  is  positive  to  the  right,  and  z  is  positive  down* 
Angles  are  measured  such  that  a  roll  to  the  right,  a  pitch 
up  and  a  yaw  right  are  positive.  Note,  however,  that 
translations  in  the  y  direction  and  roll  and  yaw  are  not 
included  in  this  model. 

Draft,  the  distance  of  the  keel  below  the  calm  water 
line,  is  the  distance  of  the  C.G.  above  the  calm  water  line, 
z,  plus  the  distance  of  the  keel  below  the  C.G.,  Zs,  or: 

Ld  =z+Zs 


(II-1) 


SIDEHALL  GEOMETRY 
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The  constant  sidewall  area  as  used  by  Gerba  and  Thaler 
yields  a  good  approximation  to  craft  motion  in  the  heave  and 
pressure  modes.  However  ,it  is  a  source  of  strong  coupling 
between  draft  and  pitch,  and  a  more  detailed  modeling  is 
necessary  to  achieve  good  agreement  with  the  Loads  and 
Motion  Program,  as  well  as  known  craft  behavior  in  pitch. 

Fig  1  shows  the  cross-sectional  detail  of  assumed 
sidewall  geometry.  The  fore  and  aft  sections  are  identical 
with  the  exception  of  deadrise  angle,  that  angle  being  60° 
at  the  bow  and  80°  at  the  stern.  The  actual  sidewall  width 
used  for  computation  of  buoyancy  forces  is  the  effective 
sidewall  width;  this  is  the  average  value  of  the  waterline 
sidewall  width  and  the  keel-line  width, Ws10  or  Ws20  (fore  or 
aft  section) .  This  is  one-half  the  draft  averaged  from  the 
C.G.  to  the  bow  and  stern  respectively,  divided  by  the 
tangent  of  the  appropriate  deadrise  angle  plus  the  keel-line 
width  Ws10  or  Ws20.  This  value  is  constrained  between  zero 
and  ws,  the  maximum  sidewall  width.   Thus: 

Ldbarl       =  Ld- (L/2-Xcg)  *Theta 

(II-2) 

Ws1  =Ldbar  1/(2*tan  (Dr  1)  ) +Ws10 

IF  (Ws1  .GT. Ws)  Ws1  =  Ws 
IF  (Ws1 .LT.O. 0) Ws1=0.0 

(II-3) 

and  for  the  stern  section: 
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Ldbar2       =Ld+ (L/2+Xcg)  *Theta 


(II-4) 


Ws2  =  Ldbar2/(2*tan (dr2)  )+Ws20 

IF  (Ws2.GT.Hs) Ws2=Hs 
IF  («s2.LT.O)  Ws2=0.0 


(H-5) 


D.   SYSTEM  DYNAMICS 

The  equations  of  the  system  dynamics  are  developed  in 
this  section  and  are  an  extension  of  the  model  of  ref  £ 5j , 
modified  as  necessary;  the  remaining  eguations  are  derived 
ficm  first  principles. 

1 •   Plenum  Pressure  and  Mass  Flow  Rate 

Plenum   pressure   and  plenum  air  mass  are  calculated 

using  the  adiabatic  law.   Using  Ld  as  the  draft  ?z    the   bow, 

1 

and   Ld  as  the  draft  at  the  stern,  the  average  draft  overall 
2 

becomes: 


Ldbar        =  (Ld+Ld)/2 

1   2 


=  (Ld- (L/2-Xcg) *tan  (Theta) 
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+  Ld+  (L/2  +  Xcg)  *tan  (Theta)  )  /2 
=Ld*Xcg*tan (Theta) 


(II-6) 


which,  when  multiplied  by  the  Plenum  area, yields  the 
submerged  Plenum  volume.  Subtracting  this  from  Plenum 
nominal  vclume  yields  the  net  Plenum  volume,  or: 

Vb  =Vn-Ab*  (Ld+Xcg*tan  (theta) ) 

(H-7) 

The  adiaratic  relationship  between  pressure  and  volume 
yields: 

Gamma 
Pb  =  Pa*  (Mb/(Vb*Rhoa)  ) 

(II-8) 

and 
Pbbar         =Pb-Pa 

(H-9) 

Vclume  flow  into  the  Plenum  is,  assuming  a  constant  fan  map: 
Qin  =N*  (QJO-Pbbar) 

(11-10) 

and  flow  cut  is: 

1/2 
Qout         =Cn*Al*  (2*Pbbar/Rhoa) 

(H-11) 
thus  the  net  mass  flow  rate  within  the  plenum  is: 
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d/    (Mb)    =Rhoa*  (Qin-Qout) 
/dt 


(11-12) 


2.   Sum  of  Forces 


Fig  2  shows  the  forces  and  their  respective  monrent 
arms  acting  on  the  XR-3.  The  sum  of  vertical  forces  yields 
the  vertical  C.G.    acceleration: 


2 

d/  2  (z)       - (Hpres+Hbuoy+Hseal+Hdamp+Hplan+W) /H 
/dt 


(11-13) 

The    individual    force   contributions   may   be   calculated 
separately  and  their  total  contributions  summed. 

a.   Pressure  Forces 

The  plenum  pressure  acting  on  the   roof  of.      the 
plenum  produces  a  pressure  force  acting  at  the  C.P.: 

Hpres        =-Ab*£>bbar 

(11-14) 


fc.   Buoyant  Forces 

For  purposes  of  calculating  buoyant  moments  in 
the  subsequent  calculations,  it  is  convenient  to  dissociate 
the  buoyant  force  into  two  components, Hbuoy 1  and  Hbuoy2,  the 
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buoyant  forces  generated  by  the  submerged  volume  fore  and 
aft  of  the  C.G.  Buoyant  force  is  proportional  to  submerged 
volume  multiplied  by  the  density  of  the  fluid  and  the 
gravitational  constant.  This  term  must  be  doubled  to  account 
for  both  sidewalls..  The  average  draft  forward  of  the  C.G. 
is: 

Ldbarl       =(Ld+Ld1/2 

-Id- (L/2-Xcg)  *tan (Theta)/2 

(H-15) 
Thus  the  forward  buoyant  force  is: 

Hfcuoyl       =-2*Rho*G*length*draft*sidewall  width 
=-2*Rho*G* (L/2-Xcg) *Ldbar1*Ws1 


(11-16) 


The  aft  average  draft  is: 

Ldbar2       =Ld+ (L/2+Xcg) *tan  (Iheta) /2 


(11-17) 


and    the    aft    buoyant   force    is: 

Hbuoy2       =-2*Rho*G* (L/2+Xcg) *Ldbar2*Ws2 


(11-13) 


and : 

Hbuoy        =Hbuoy 1 *Hbuoy2 


(11-19) 


c.   Eow  and  Sternseal  Forces 
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Fig  3  shows  the  detail  of  the  assumed  bowseal  and  sternseal 
geometry  and  water  surface  interactions.  The  seal,  is 
assumed  to  trail  from  the  upper  hinge  point  at  a  constant 
angle  of  31.6°  until  it  contacts  the  waterline.  From  that 
point  on  it  is  assumed  to  lie  tangent  to  the  waterline.  The 
length  of  bowseal  in  contact  with  the  water  may  be  seen  to 
be: 

Xseall       =Ld1/cos(31.60) 

=  1.  65*  (Ld-  (L/2-Xcg)  *tan  (Theta)  ) 

(11-20) 

Pbbar  exerts  a  force  upon  this  wetted  area  which  is  resisted 
by  an  upward  reaction  from  the  water  surface.  This  reaction 
is  transmitted  to  the  craft  as: 

Hseall       =-Pbbar*Width*Xseal1 

(11-21) 

As  may  be  seen  from  Fig  3,  the  sternseal 
gecmetry  is  basically  identical  to  that  of  the  bowseal; 
however,  the  mechanism  for  force  generation  is  different  in 
that  the  sternseal  rides  atop  an  exhausting  jet  of  plenum 
air.  Therefore,  the  pressure  differential  is  not  that  of 
Plenum  pressure  over  atmospheric,  but  over  the  reduced 
pressure  of  the  high-velocity  jet.  This  is  a  complicated 
function  of  leakage  area  and  other  factors,  but  aay  be 
assumed   approximately   constant.     This    model    assumes 

2 

approximately   2.0   lbf/ft  differential  pressure  between  the 

sternseal  and  plenum  pressure. 

The  length  of  sternseal  in  contact  with  the 
water  may  be  seen  to  be: 
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Xseal2       =Ld2/cos (3 1 . 6°) 

=  1 . 65*  (Ld+ (L/2+Xcg) *tan (Theta) ) 

(11-22) 

the  sternseal  force  is  then  the  differential  pressure, Sealp, 
times  the  sternseal  area  ,or: 

Hseal2       =-Sealp*Width*Xseal2 

(H-23) 

d.   Damping  Forces 

Damping  forces  are  developed  for  information 
only,  as  Taylor's  series  expansion  of  these  forces  about  the 
zero-rate  equilibrium  point  causes  these  equations  to  vanish 
in  the  linear  model. 

A  body  at  motion  within  a  fluid  experiences  a 
frictional  force: 

F  =-C*  (Rho/2)  *V*|  V|*S 

f  f 

where  C  is  a  non-dimensional  friction  coefficient  based  upon 
f 

shape,  roughness,  and  Reynold's  Number.   Consistent  with  the 

six-degree-of-freedom  Loads  and  Motion  program,   C  is   taken 

to  be  .9.   The  vertical  velocity  is  taken  to  be  the  vertical 

velocity  of  the  C.G..,   Zdot,  plus  the  tangential  velocity  of 
the  angular  velocity,  which  varies  along  the  hull  as: 
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V  (x)  =Zdot>  (Xcg-x)  *Omega 


(11-24) 


The  differential  surface  wetted  area  also  is  a  function  of 
x,  as  the  depth  of  immersion  varies  along  the  hull,  so  that: 

dS  =Ws*(Ld-(x-Xcg)*tan  (Theta)  ) 

(H-25) 

When  multiplied  by  the  absolute  value  of  velocity,  times 
itself,  times  the  differential  wetted  area,  the  resulting 
expression  may  be  integrated  numerically  along  the  null. 
However,  this  is  most  wasteful  of  CPU  time,  and  the  absolute 
value  term  renders  this  expression  intractable  to  closed 
analytical  solution..  In  addition,  the  term  is  quite  small, 
and  in  fact  must  be,  if  pitch  angles  are  to  remain  within 
the  small  angle  assumptions,  and  no  sidewall  gapping  occur. 
Therefore,  this  term  is  approximated  by  computing  the 
average  velocity  fore  and  aft  of  the  center  of  gravity,  and 
the  average  wetted  surface  fore  and  aft,  previously  computed 
as  Ws1*Ldbar1  and  Ws2*Ldbar2.  For  purposes  of  pitch 
calculations  to  be  discussed  subsequently,  these  two  forces 
are  assumed  to  act  at  two-thirds  the  length  fore  and  aft  the 
C.G.;  this  yields  closer  agreement  to  the  numerical 
solution  by  weighting  the  higher  velocities  which  occur  at 
the  bow  and  stern.   Vbarl  and  Vbar2  may  now  be  defined: 

Vbarl        =Zdot-  (L/2-Xcg) *Omega/2 

(11-26) 

Vbar2         =Zdot+  (L/2 +Xcg) *0mega/2 

(H-27) 
Since   there  are  fcur  surfaces  forward  and  aft  of  the  center 
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of  gravity  the  damping  terms  evaluate  as: 

Hdamp1=-2C*Rho*Vbar1*| Vbar 1|*Ldbar1*  (L/2-Xcg) *Ws1 
f 


(11-28) 


Hdamp2=-2C*Rho*Vbar2* |Vbar2 | *Ldbar2* (L/2+Xcg) *Ws2 


(11-29) 


e.   Planing  Forces 

Initial  time  domain  validation  studies  conducted 
with  the  model  developed  from  the  eguations  of  this  chapter, 
and  the  Loads  and  Motions  program,  disclosed  a  serious 
discrepancy  between  the  two  models  in  natural  pitch 
freguency.  Specifically,  the  Loads  and  Motions  program 
predicted  a  natural  pitch  freguency  of  3.2  rad/sec  at  10 
kts,  and  5.4  rad/sec  at  20  kts,  compared  with  a  constant 
pitch  freguency  of  3'.  3  rad/sec  for  the  linear  model.  This 
was  taken  to  indicate  the  presence  of  a  velocity-dependent 
source  of  theta  feedback.  Examination  of  ref  [2.J  disclosed 
that  Leo  and  Boncal  had  indeed  utilized  an  added  mass  effect 
which  had  the  effect  of  a  planing  force,  in  their 
modifications  of  the  Loads  and  Motions  program. 

A  body  in  motion  in  a  fluid  experiences  a 
lifting  force  proportional  to  the  sguare  of  the  velocity  and 
directly  proportional  to  the  angle  of  attack,  or: 

2 
Hplan        =-C*Theta*Rho*V*S/2 

la 

(11-30) 
where  S  is  the  surface  area  of  the  body,  or: 
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=  (Ws1*Ws2)  *L/2 


(11-31) 

and  C    is  the  lift-curve  slope.   This  term  is  not   directly 
la 

available,   and   is   a   complex   function   of   the   sidawall 

gecmetry.  A  value  of  1.5  gave  good  agreement  with  the 
natural  frequency  of  the  Loads  and  Motions  program,  but  a 
less  empirical  approach  to  this  value  is  recommended  for 
future  studies. 


The  planing  force  is  assumed  to  act  at   the   aft 
quarter-length  point. 

3 .   Sum  of  Moments 

The  forces  described  above  are  multiplied  by  their 
respective  mcment  arms  as  shown  in  Fig  2  to  yield  the 
pitching  moments  about  the  center  of  gravity.  The  angular 
acceleration  about  the  C.G.  is  then: 

Alfa  = (Ppres+Pbuoy+Pseal+Pdamp+Pplan) /Iyy 

(11-32) 

a.   Pressure  Moment 

The  pressure  exerts  a  moment  about  tne  center  of 
gravity,  equal  to: 
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Ppres        =-Hpres*  (Xcp-Xcg) 

(11-33) 

b.   Buoyant  Moments 

The  buoyant  forces  due  to  the  immersed  volumes 
fore  and  aft  of  the  center  of  gravity  exert  a  pitching 
mcment  alcng  their  moment  arms,  which  are  approximately  (for 
small  angles)  the  centroid  of  the  rectangle  from  the  center 
of  gravity  to  the  bo.w  and  stern: 


Pbuoyl       =-.5* (L/2-Xcg) *Hbuoy1 


(11-34) 


Pbuoy2       =.5*(L/2+Xcg) *Hbuoy2 


(11-35) 


c.   Bow  and  Sternseal  Moments 

The  bowseal  acts  along  a  moment  arm  that  is 
formed  by  the  distance  from  the  C.G.  to  the  centroid  of  the 
wetted  area  Aseal,  or: 

Pseall       =-Hseal1*(L/2-Xcg-Xseal1/2) 

(11-36) 

The  sternseal  force  acts  on  a  moment  arm  that  is 
formed  by  the  distance  from  the  C.G.  to  the  csatroid  of  the 
sternseal  wetted  area,  or: 


U1 


Pseal2  =Hseal2*(L/2+Xcg+Xseal/2) 

(11-37) 

d.  Damping   Moments 

The  damping  moments  in  pitch  are  the  damping 
forces  in  heave  multiplied  by  their  respective  moment  arms, 
two-thirds  the  distance  from  the  C.G.  to  the  bow  and  stern 
respectively : 

Pdampl       =-  (2/3) *(L/2-Xcg)  *Hdamp1 

(11-38) 

Pdamp2  =  (2/3)  *  (L/2  +  Xcg)  *Hdamp2 

(11-39) 

e.  Planing   Moments 

The  planing  moment  is  the  planing  force  acting 
at   the   aft    quarter-length    point,    or: 

Pplan  =  (L/4  +  Xcg)  *Hplan 

(11-40) 


E.   INITIAL  CONDITIONS 


As  in  the  Loads  and  Motion  program,   initial   conditions 
must  be  carefully  selected  to  ensure  steady  state  conditions 
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at  the  beginning  of  a  run.  Since  speed  does  not  enter  into 
this  program,  one  of  the  major  problems  of  the  Loads  and 
Motions  program  is  obviated.  However,  the  craft  in  this 
model  is  particularly  sensitive  to  pitch,  and  once  excited 
in  this  mode,  is  ver.y  lightly  damped. 

Selection  of  initial  conditions  proceeds  as  follows: 

1.  pitch  rate  and  vertical  velocity  are  defined  as  zero; 

2.  a  value  of  Pbbar  is  selected  such  that  plenum  air  mass 
flow  rate  is  set  to  zero; 

3.  a  value  of  draft  is  calculated  such  that  weight  and 
pressure  forces  are  balanced  by  buoyant  and  planing 
forces; 

4.  a  value  of  Theta  is  calculated  such  that  pressure 
moments  are  balanced  by  buoyant  and  planing  moments; 

5.  given  the  Plenum  volume  defined  by  the  equilibrium 
values  of  draft  and  Theta,  and  the  equilibrium  value  of 
Pbbar,  the  Plenum  air  mass  may  be  calculated 
adiatatically . 


More  detail  on  the  exact  procedures  for  the  selection  of 
equilibrium  conditions  may  be  found  in  subsequent  chapters 
dealing  with  the  computer  simulation  of  this  model. 
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III.   LINEAR  HEAVE-PITCH  MODEL 


The  non-linear  equations  of  motion  developed  in  the 
preceding  chapter  form  the  basis  of  a  linearized  heave-pitch 
modelr  valid  for  disturbances  about  its  equilibrium  point, 
defined  as  that  point  aoout  which  all  rates  are  zero.  The 
linear  model  permits  the  use  of  a  variety  of  powerful 
techniques  for  parametric  analysis,  such  as  root  locus;  it 
furthermore  simplifies  the  design  of  an  automatic  control 
system,  which  could  serve  to  further  limit  excursions  about 
the  operating  point. 


A.   TAYLOR'S  SERIES  EXPANSION 


For   small    excursions 
non-linear  relationships: 


about    an   operating   point, the 


f  =F(x,y) 

g  =G(x,y) 

are  approximately: 


f  +  df 
0 


=  F  (x,y)  +(dF/dx)  *Dx  +  (dF/dy)  *Dv 
0  0 


go   d9 


=G(x  ,y  ) + (dG/dx) *Dx+ (dG/dy) *Dy 
0   0 


where  the  derivatives  of  F  and  G  with  respect  to  x  and  y  are 

evaluated   at   the  operating  point,  x  and  y.   Cancelling  the 

0      0 

equilibrium    point    from    both    sides    and      dividing      through      by 
time   yields: 
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a/   (f) 
/at 


=  {dr/dx)*  d/    (x)  +  (d?/dy)  *  d/    (y) 
/at  /dt 


d/    (q) 

/dt 


=  (dG/dx)*  d/    (x)  +  (dG/dy)*  d/    (yj 

/dt  /It 


Applying  these  principles  to  the  non- linear  model  of  the 
preceding  chapter,  the  steady  stare  operating  point  ie 
defined  as  before  for  the  equilibrium  values  of  the  state 
variables. 


3.   PL£NUM  GAUGE  ?R£33U?.i:  DEHIVATIVES 


Plenui  gauge  pressure  has  jaen  ueen  defined  as 


CDar 


=Pb-?a 


(II-9) 


and  since  Pa  is  a  constant 


dPtbar 


=  dPo 


Pb,    the    plenum    total    pressure,    udo    been    defined   as 


Pd  (0) 


=  ?i*  (:•]  j  (0)  /  (RhOa-*lD  (J)  ) 


baifima 


Therefore  expanding  arcund  the  equilibrium  values   of 
and  Hb  (0)  : 


(II-6J 
Vn  (0) 


Eb  (0)  +  dPb     =Pa*  (Mb  (0)  /  (3hoa*Vb  (0)  ) 


Gaana 
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+  (Gamma/Mb  (0)  )  *Pa*(Mb  (0)  /  (Bhoa*Vb  (0)  )  *dM£ 

Gamma 
- (Gamma/Vb (0) ) *Pa* (Mb (0) / (Rhoa*Vb (0) ) *dVb 


(III-1) 

or,  after  substituting  the  definition  of  Pb(0)  and 
cancelling  the  equilibrium  values  from  both  sides: 

dPb  =Gamma*  (Pb(0) /Kb(0)  )*dab 

-Gamma*  (Pb  (0)  /Vb  (0)  )  *dVb 

(HI-2) 

Vt  (0)  may  te  further  expanded, since  it  has  been  defined  in 
the  preceding  chapter  as  the  empty  plenum  volume  less  the 
net   submerged   volume,,    following    the   definition   of    ref    jjSJ : 

Vb(0)  =Vn-Ab*  (z(0)  +Zs+Xcg*Theta  (0)  ) 

(II-7) 

and 
dVb  =-Ab*dz-Ab*Xcg*dTheta 

(III-3) 

Substituting  the  expansion  of  incremental  volume,  Equation 
(III- 3) ,  into  the  expansion  for  incremental  pressure, 
Equation  (III-2) ,  and  dividing  by  differential  time  yields: 


d/    (Pbbar)=Gamma*  (Pb  (0) /Mb(0)  )  *  d/    (Mb) 
/dt  /dt 
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+Gainaa*Ab*  (?b  (J)/Vb  (0)  )  *    a/         (z) 

/at 

+  Gdiama*Ab*Xcg*  (pb  (0) /Vb(0)  )  *    a/       (Theta) 

/dt 


For    mors    compact    notation    the 
coefficients  are  defined  such  that: 


(111-4) 

following   influence 


Dpbz=Gamma*Ab* (Pb (0)  /Vo (0)  )  ; 


Dpbth=Gaaaa*Axj  *Xcg*  (Pb  (0)  /Vb  (0)  )  ; 


Dpbmb=Gamma*  (Po  (•))  /Mb  (0)  )  ; 


and    the    pressure    iifferential   equation    becomes; 


d/  (Pttar)=0?cz*    d/         (z) 

/dt  /dt 


+Dpbth*    d/  (Theta) 

/dt 


+  0pb:nc*    1/  (.ID) 

/it 


C.       PLENUM    AIR    MASS    FLC«     RATE    DERIVATIVES 


(I1I-5) 


(III-6) 


(Tir-7) 


(III-8) 
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The  plenum  air  mass  flow  rate  has  been   defined   in   the 
preceding  chapter  as; 


d/ 

/dt 


(Mb(0))  =fihoa*  (Qin  (0)  -Qout  (0)  ) 


and 


d(  d/    (Mb)  )=Rhca*dQin-Rhaa*dQout 
/dt 


Qin  has  been  defined  as: 

Qin(0)        =N*  (Q^O-Pbbar  (0) ) 


(11-12) 


(III-9) 


so  that 
dQin 


=-N*dPbfcar 


(11-10) 


and  Qout  has  been  defined  as: 


Gout  (0) 


=Cn*Al*  (2*Pbbar  (0)  /Ehoa) 


1/2 


(111-10) 


and 


dQout 


(11-11) 


-1/2 

(Cn*Al/Rhoa)  *  (2*Pbbar  (0)  /Rhoa)  *dPbbar 


(111-11) 

Combining   the   two   expressions   for  incremental  volumetric 
flew  into  and  out  of  the  plenum,  Equations  (111-10,11)  ,  with 
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the   expression  for  .for  incremental  mass  flow  rate,  Equation 
(III-9)  ,  and  dividing  by  differential  time  yields: 

2  -1/2 

d/  2  (Mb) =- (N*Rhoa+Cn*Al* (2*Pbbar (0) /Rhoa) ) *  d/    (Pbtar) 
/dt  /dt 

(111-12) 

The  influence  coefficients  may  be   defined   by   substituting 
Equation  (III-8)  into  Equation (111-12)  and  eguating  terms: 

-V2 

Dmbz         =- (N*Rhoa+Cn*Al*  (2*Pbbar  (0) /Rhoa) ) *dpbz 


(111-13) 

-1/2 
Dmbth  =-(N*Rhoa+Cn*Al*(2*Pbbar (0) /Rhoa) ) *dpbth 


(111-14) 

-1/2 

Dmbmb  =- (N*Rhoa+Cn*Al* (2*Pbbar  (0) /Rhoa) ) *dpbmt 


(111-15) 

Thus   the   linearized   differential  equation  for  plenum  mass 
flow  rate  may  be  written: 

2 
d/  2  (Mb)    =Dmbz*  d/    (z) 
/dt  /dt 


♦  Dmbth*  d/    (Theta) 
/dt 


♦Dmbmb*  d/    (Mb) 
/dt 


(111-16) 
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D.       HEAVE    DERIVATIVES 


The    heave   equation   has   been   defined   as: 

2 
d/   2     (z)  =  (Hpres+Hbuoy+Hseal+Hdamp+Hplan+tf)  /M 

/dt 


(11-13) 

Damping  forces  are  neglected  in  the  linear  model. We  shall 
consider  the  various  components  of  heave  force  separately, 
and  superimpose  the  components  upon  completion  of  all 
derivations. 


1  •   M.^§v  e  Force  I£ ue  To  Pressure 

The  heave  force  due  to  pressure  has  been  defined  as: 
Hpres        =  -Ab*P.bbar 


(11-14) 


By  inspection,  its  influence  coefficients  become: 

Dhpz=-Ab*Dpbz 


(111-17) 


Dhpth=-Ab*Dpbth 


(III-18J 


Dhpmb=-Ab*Dpbmb 


(111-19) 
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2 •   H£a vg  Forces  Due  To  Buoyancy 

The  buoyant  heave  force  equation  from  the 
chapter  has  been  defined  as: 

Hbuoy        =Hbuoy 1+Hbuoy2 
=-2*Rho*G*Ws1*(L/2-Xcg)*  (Ld-  (L/2-Xcg)  *Theta/2) 
-2*Rho*G*Ws2*  (L/2  +  Xc^)  *  (Ld+  (L/2+Xcg)  *Tneta/2) 


preceding 


(11-15, 16,  17,18) 


which  may  be  rearranged: 
Hbuoy= 

-2*Rho*G*  <Ws  1*  (L/2-Xcg)*ws2*  (L/2+Xcg)  )  *  (z  (0)  +Zs) 


2  2 

+2*Rho*G*  (Ws  1*  (L/2-Xcg) -Ws2*  (L/2+Xcg) ) *Thexa/2 


and    the    influence   coefficients   may    be    written    by   inspection: 
Dhbz  =-2*Rho*G*  (Ws1*  (L/2-Xcg) + Ws2*  (L/2+Xcg) ) 

(111-20) 


Chbth 


2  2 

=Rho*G*  (Ws1*  (L/2-Xcg)  -Ws2*  (L/2  +  Xcg)  ) 


(111-21) 


3 .       H<saye    Force   due   to    Bow    and   S  ternseal 
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The   expression  for  the  bowseal  heave  force  has  been 
defined  as: 

Hseall  (OJ     =- Width*Xseal  1  (0)  *?bbar  (0) 

(H-21) 
which  for  small  perturbations  in  Xseall  and  Pbbar  becomes 
dHseall       =-Width*Pbbar  (0) *dXseal  1 

-Width*Xseal1 (0) *dPobar 

(111-22) 
The  expression  for  Xseall  has  been  defined: 
Xseall  (0)     =1.65* (z(0)  +Zs- (L/2-Xcg) *Theta) 

(11-20) 

whj.ch  for  small  perturbations  around  z(0)  and  Theta  (0)  is 
approximately 

dXseall       =1 . 65*dz-1 .65* (L/2-Xcg) *dTheta 

(111-23) 

Substituting  these  relations  along  with  the  relations  for 
small  perturbations  in  Pbbar  due  to  z, Theta  and  Mb  from 
Equations  (111-5,6,7)  yields: 

dHseall      =-  Width*  (1 .  6  5*Pbbar  (0)  +Xseal1  (0)  *Dpbz)  *dz 

-Width* (- 1 .65* (L/2-Xcg) *Pbbar  (0) 

♦Xseall (0) *Dpbth) *dTheta 
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-Width*Xseal1  (0)  *Dpbmb*dMb 

(111-24) 
The   influence   coefficients   then   become: 

Dhszl  =-»id±h*(1 .65*Pbbar  (0)  +Xseal1  (0)  *Dpbz) 

(111-25) 

Dhsthl       =-Width*(-1.6  5*  (L/2-Xcg)  *Pbbar  (0) 

+Xseal1 (0) *Dpbth) 

(111-26) 

Dhsmb         =-Width*Xseal1 (0)  *Dpbmb 

(111-27) 

The  expressipn  for  the  stern  seal  has  been  defined: 
Hseal2       =-Sealp*Width*Xseal2 (0) 

(11-23) 

Xseal2(0)      =1.65*  (Ld  (0)  +  (L/2+Xcg)  *tan  (Theta  (0)  )  ) 

(11-22) 
Thus  for  small  perturbations  in  draft  and  Theta: 
dHseal2      =- 1 .65*sealp* Width*dz 

-1.65*(L/2+Xcg) *Sealp*Widt h*dTheta 

(111-23) 
Thus  the  steinseal  influence  coefficients  become: 
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Ehsz2        =-1 .65*Seal?*Width 

(111-29) 
Ehsth2       =-1 .65*  (L/2  +  Xcg) *Sealp*Width 

(111-30) 

**  •      J§^£i  l£iL£$  4ue  to    Planing 

The  planing  force  is  linear  in  Theta  ,  or: 

2 
Hplan        =-C  *Sho*V*(Ws1+Ws2)*L*Theta 

la 

(11-30,31) 
and  the  influence  coefficient  may  be  written  by  inspection: 

2 
Dhplth        =-C*Rho*V*(Ws1+Ws2)  *L 

la 

(111-31) 
5.   Sum  Cf  Forces 


The  forces  previously  derived  may  be  super  imposed  so 
that: 


2  2 

d/  2  (z(0))+d(  d/  2  (z))  = 

/dt  /dt 
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(Hpres  (0)  +  Hbuoy  (0)  +Hseal(0)  ♦Hplan(O)  )/M 
+  (Dhpz+Dhbz  +  Dhsz)  *dz/M 
+ (Dhpth+Dhbth+Dhsth+Dhplth) *dTheta/M 
+  (Dhpmb+Dhsmb) *dflb/H 

(111-32) 
The   influence   coefficients    may    now    be   defined: 

Dzz=  (Dhpz+Dhbz+Dhsz)  /M 

(111-33) 
Dzti=  (Dhpth+Dhbth+Dhsth+Dhplth)  /M 

(111-34) 
Dzmb=  (Dhpmb  +  Dhsmb) /M 

(111-35) 

Hemoving  the  equilibrium  values  and  dividing  through  by 
differential  time  yield  the  differential  eguation  of  motion 
in  the  z-direction: 

3 
d/  3  (z)     =Dzz*  d/    (z) 
/dt  /dt 


+Dztn*  d/    (Theta) 
/dt 


+Dzmb*  d/    (Mb) 
/dt 


(111-36) 


E.   PITCHING  MOMENT  DEBIVATIVES 
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The    pitching    moment  equation    is 

2 

d/  2  (Tbeta)  = (Ppres+Pbuoy  +Pseal+Pdamp+Pplan) /Iy y 
/dt 


(11-32) 

As  in  the  previous  case  of  heave  forces,  damping  moments  are 
neglected  in  the  linear  model;  the  various  moments  are 
considered  separately  and  superimposed  after  derivation. 

1 •   li^ch  Moment  due  to  Pressure 

The  pressure  moment  has  been  defined  as 
Ppres         =  (Xcg-Xcp)  *Hpres 

(11-33) 
The  influence  coefficients  may  be  written  by  inspection: 

Dppz=  (Xcg-Xcp)  *Dhpz 


(111-37) 


Dppth= (Xcg-Xcp) *Dhpth 


(111-38) 


Dppmb=  (Xcg-Xcg)  *Dhpmb 


(111-39) 


2«   Eu°X§.D.t  Moments 

The  buoyant  (moment  has  been  defined 
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Pbuoy        =Pbuoy 1+Pbuoy2 

=- (L/2-Xcg) *Hbuoy1/2* (L/2  +  Xcg) *Hbuoy2/2 

2 
=Rho*g*  (Ws1*  (L/2-Xcg) * (Id- (L/2-Xcg)  *£heta/2) 

2 
-Ws2*  (L/2  +  Xcg)  *  (Ld+  (Ld+  (L/2  +  Xcg)  *Theta/2)  )  ) 

(11-34,35) 

and  the  influence  coefficients  may  be  defined: 

2  2 

Dpbyz        =Rho*G*  (Ws1* (L/2-Xcg)  -Ws2* (L/2  +  Xcg)  ) 

(111-40) 


3  3 

Dpbyth         =-Rho*G*(Ws1*  (L/2-Xcg) +  Ws2*  (L/2+Xcg) ) /2 


(III -41) 


3.   Sow  and  5  ternseal  Moments 

The  towseal  contributes  a  moment 
Pseall         =-Hseal 1* (L/2-Xseal1/2-Xcg) 


(11-36) 


which  for  small  disturbances  ,  is  approximately 
dPseall        =- (L/2-Xseal1  (0) /2-Xcg) *dHseal1 
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♦Hseall (0) *dXseal1/2 


(111-42) 


Substituting   the   definition  cf  dXseall,  Equation  (111-23)  , 
and   the   influence   coefficients    of    Hseall,    Equations 
(111-25,26, 21)    yields: 

dEseall       =  <-  (L/2-Xseal  1  (0)  /2-Xcg)  *Dhsz1 

+1.6  5*Hseal1  (0)  /2)  *dz 

♦  (-  (L/2-Xseal1  (0)/2-Xcg)  *Dhsth1 

-  1.65  (L/2-Xcg)  *Hseal1  (0)  /2)  *dTheta 

-(L/2-Xseal1  (0)  /2-Xcg)  *Dhsmt*dMb 


(III-43J 


The  influence  coefficients  may  now  be  defined: 

Dpszl        =- (L/2-Xseal1  (0) /2-Xcg) *Dhsz1 

+1.65*Xseal1  (0)  /2 


(111-44) 


Dpsthl  =   -  (L/2-Xseal1  (0) /2-Xcg) *Dhsth1 

-1.65*(L/2-Xcg)  * Hseall  (0)  /2 


(111-45) 


Epsmb  =- (L/2-Xseal1  (0) /2-Xcg) *Dnsmb 


(111-46) 


The  sternseal  contributes  a  moment: 
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Pseal2  (0)     =Hseal2* (L/2+Xseal2 (0)  /2+Xcg) 


which  for  small  disturbances  is  approximately 
dPseal2      = (L/2+Xseal2 (0) /2+Xcg) *dHseal2 

+Hseal2 (0) *dXseal2/2 


(11-37) 


and   substitution   from   the    sternseal    force 
coefficients,  Equations  (111-29,30),  yields: 


dPseal2 


(111-47) 

influence 


=  ( (L/2  +  Xseal2 (0) /2  +  Xcg)  *Dhsz2 
+  1.65*Hseal2(0)  /2)  *dz 

+  ( (L/2+Xseal2 (0) /2  +  Xcg)  *Dhsth2 
+  1.65*(L/2  +  Xcg) *Hseal2 (0) /2) *dTheta 


(III-48J 

and  the  stern   seal   pitch   influence   coefficients   may   be 
defined: 

Dpsz2=(L/2+Xseal2(0) /2+Xcg) *Dhsz2 

+1.65*Hseal2  (0)  /2 


(111-49) 


Dpsth2=(L/2  +  Xseal2 (0)  /2  +  Xcg) *Dhsth2 
+  1.65* (L/2+Xcg) *Hseai2 (0)  /2 


(111-50) 


*•  •   Pitch  Moment  due  to  Plani n^ 
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The  planing  moment  is  linear  in  Theta: 
Pplan        =  (L/4+Xcg)  *Hplan 

(11-40) 
and  the  influence  coefficient  is: 
Dpplth       =  (L/4*Xcg)  *Dhplth 

(111-51) 

5 .   Sum  C J  Moments 

The   individual   moment  contributions   may   new   be 
summed: 

Alfa  (0)  +dAifa-  (Ppres  (0)+Pbuoy  (0)  +Pseal  (0)  +Pplan  (0)  )  /Iyy 

+ (Dppz+Dpbyz+Dpsz) *dz/Iyy 

+  (Dppth+Dpbyth+Dpsth+Dpplth) *DTheta/Iyy 

+  (Dppmb+Dpsmb) *dMb/Iyy 

(111-52) 

The  influence  coefficients  may  now  be  defined 

Dp z= (Dppz+Dpbyz+Dpsmb) /Iyy 

(111-53) 
Dpth=  (Dppth  +  Dpbyth  +  Dpsth  +  Dppith) /Iyy 

(1 11-54) 
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Dpmi)=  (Dppmb+Dpsmb)/Iyy 


(111-55) 


and  the  differential  equations  of  motion  for  angular  rare  of 
acceleration  may  be  written,  after  dividing  through  by 
differential  time  and  cancelling  the  equilibrium  values: 

3 
d/  3  (Theta)=Dpz*  d/    (z) 
/dt  /dt 


+Dpth*  d/    (Theta) 
/dt 


♦Dpmb*  d/    (Mb) 
/dt 


(111-56) 


F.   STATE  SPACE  REPRESENT ATION 


Fig  4  is  the  linear  model  as  derived   in   this   chapter 
presented  in  matrix  form  as: 

x  =  A*x 

Examination  of  the  eigenvalues  of  the  9-by-9  matrix  shows 
four  eigenvalues  at  the  origin  (0±0j)  which  represent  the 
integrated  output  equations  for  zr  Theta,  Mb,  and  Pbbar;  the 
remaining  eigenvalues  consist  of  a  complex  pair,  located  in 
the  vicinity  of  -30±30j,  and  a  real  eigenvalue  in  the 
vicinity  of  -.5,  as  j.n  ref  ],5J,  and  an  additional  complex 
pair  near  0±4j,  representing  the  very  lightly-dammed  pitch 
response . 
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Figure  4  -   STATE  SPACE  REPRESENTATION  OF  THE  LINEAR  SYSTEM 
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IV.   TIME  DOMAIN  VALIDATION 


The  linear  heave-pitch  XR-3  model  was  validated  against 
the  modified  Oceanics  XR-3  Loads  and  Motions  Program  of  ref 
£3  J  prior  to  investigations  in  the  frequency  domain.  The 
linearized  program,  which  appears  at  the  end  of  the  text, 
was  simulated  using  the  Continuous  System  Modeling 
Program,  version  III/.  This  is  basically  a  FORTRAN  language 
with  extremely  powerful  variable-step  integrators  and 
versatile  output  structure.  Readers  familiar  with  FORTRAN 
should  find  little  difficulty  following  the  CSMP  III 
programs  developed  for  this  simulation;  for  details  of  this 
simulation  language,  see  ref  f6J . 

The  Time  Domain  Program  consists  of  three  sections; 

1.  Initial  Conditions;  this  subsection  identifies  the 
equilibrium  operating  point  for  a  given  set  of  craft 
physical  parameters. 

2.  Linearization;  the  partial  derivatives  of  the 
non-linear  state  equations  derived  in  the  preceding 
chapter  are  evaluated  at  this  particular  operating 
point. 

3.  Dynamics;  the  weight  and  C.G.  acceleration  are 
perturbed  at  time  zero  by  a  step  weight  change,  and  the 
resultant  linear  state  trajectories  in  time  are 
computed.  Note  that  the  linear  model  is  not 
self-starting;  an  initial  C.G.  acceleration  must  also 
be  applied. 


66 


A.   INITIAL  CONDITIONS 

Initial  conditions  involve  the  solution  of  the 
non-linear  system  equations  of  Chapter  II  for  the 
equilibrium  operating  point  defined  at  the  end  of  the 
chapter. 

1 •   ^i^num  Pressure 

A  value  of  plenum  pressure  is  selected  so  that 
plenum  mass  flow  rate  is  set  to  zero,  or,  from 
Equations  (11-10,11)  : 

gin  =Qout 

1/2 
CiO-Pbbar  = (Cn*Al/N) *  <2*Pbbar/Rhoa} 


(IV-1) 


which  is  a  quadratic  in  Pbbar: 

2  2  2 

Pbbar-  (2*QiC+  (Cn* Al/N)/Rhoa) *Pbbar+QiO=0 


(IV-2) 

The  larger  solution  of  this  quadratic,  corresponding  to  the 
physically  unrealizable  condition  of  reverse  flow  through 
the  fans,  is  rejected. 

2-   ££§ft  and  £itch  Angle 

The  solution  for  equiliorium  values  of  draft  and 
pitch,  given  the  previously-computed  value  of  Pbbar, 
requires    an    iterative   technique,   due   to   their   close 


(1 


coupling.  An  initial  value  of  draft  is  first  computed 
assuming  2ero  pitch  and  no  bow  or  stern  seal  contribution. 
The  ncn- linear  state  eguations  are  then  solved  for  residual 
forces  and  moments  due  to  all  sources.  A  very  small 
perturbation  in  draft  and  pitch  angle  of  one  percent  with  a 
minimum  default  magnitude  of  .001  is  applied  and  tae 
differential  change  in  residual  forces  and  moments  with 
respect  to  draft  and  Theta  are  computed.  Bairstow»s  method 
for  the  zeroes  of  functions  of  two  variables,  as  described 
in  ref  [7] ,  is  used  to  drive  the  residuals  to  zero.  This 
particular  problem  is  somewhat  slow  to  converge;  convergence 
is  enhanced  by  amplifying  the  computed  changes  in  draft  and 
pitch  by  1.9,  but  after  two  hundred  iterations  residuals  on 
the  order  cf  ten  lif  and  ten  ft-lbf  may  still  be  observed. 
These  residual  forces  and  moments  are  applied  to  the  initial 
values  of  Z2dot  and  Alfa. 

3 -   Plenum  Air  Mass 

The   eguilibrium   values   of   draft   and  pitcn  angle 
define  a  unigue  plenum  volume: 


Vb  =Vn-Ab*  (Ld+Xcg*tan  (Theta)) 


(II-7) 


which,  when  combined  with  the  equilibium  pressure,  yields   a 
unigue  adiatatic  plenum  air  mass,  from  Equation  (11-8}  : 

1/Gamma 
MbO  =Rhoa*Vb* (Pbbar/Pa) 


(IV-3] 


The    solutions    to    these   eguations   define   the 
equilibrium  values  of   Pb  (0)  ,   Ld  (0)  ,   Theta  (0),   and   iib(0) 


68 


about   which   the   partial   differential   equations   of   the 
preceding  chapter  may  be  evaluated. 


B.   LOADS  AND  MOTIONS  PROGRAM  OPERATING  CONDITIONS 


The  XR-3  Loads  and  Motion  program  was  run  at  6500  lb 
weight,  C.G.  S.825  ft.  ahead  of  the  aft  transom  at  20  kts. 
Simulation  time  was  30.0  sec.  for  the  purpose  of  obtaining 
reasonably  quiescent  values  of  pitch  angle,  draft,  and 
Plenum  pressure.  In  order  to  obtain  identical  operating 
conditions,  it  was  necessary  to  adjust  various  constants  in 
the  linear  model  until  agreement  was  reacned  with  the 
steady-state  operating  conditions  of  the  Loads  and  Motion 
program.   Specifically,  these  parameters  were: 

1.  seal  leakage  area,  hi,    adjusted  until  the  two   programs 
agreed  in  Plenum  pressure 

2.  Planing  lift-curve  slope   coefficient,   C    ,   adjusted 

la 

until    the    two   programs   agreed   in   natural   pitch 
f reguency ; 

3.  minimum   sidewall   width,   WsO,   adjusted  until  the  two 
programs  agreed  in  draft; 

4.  Residual   pitching   moment  ,Mf0,  adjusted  until  the  two 
programs  agreed  in  steady-state  pitch  angle. 

The  results  of  these  triai-and-error  iterations  yielded 
the  following  table  pi   values: 
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TABLE  OF  OPERATING  PARAMETERS 


LINEAR  SYSTEM 

LOADS  AND  MOTIONS 

PARAMETEfi 

AJFECTED 
V15IABLE 

STEADY-STATE 
VALUE 

STEADY-STATE 
VALUE 

Al=.390 

C=1.5 
la 

Ws0=.27 
Mf0=-1875 

Ptbar 

n 

F 
Ld 

Iheta 

2 
24.85  lbf/ft 

5.02   rad/sec 

8.497  in. 
1.21   deg 

2 
24.84  lbf/ft 

5.46   rad/sec 

8!.55   in. 
1.199  deg 

C.   COMPARISON  OF  RE'SULTS 

The  Loads  and  Motions  program  was  run  3.0  sees  with  the 
initial  conditions  as  obtained  above,  but  with  a  weight  of 
6400  lbs.  The  linear  model  was  simultaneously  run  with  a 
100  lb  step  weight  reduction;  data  from  the  Loads  and 
Motions  program  was  subsequently  entered  into  the  CSMP  III 
linear  program  in  tabular  form,  for  simultaneous  plotting 
and  error  computation.  For  the  purposes  of  comparison,  error 
was  defined  as  the  difference  between  the  data  from  the  two 
programs,  divided  by  the  average  value  of  the  two.  This 
procedure  significantly  reduced  the  probability  of 
inadvertent  division  by  zero,  as  compared  with  normalizing 
the  difference  with  respect  to  one  system  or  the  other. 

1  •   Planum  Pressure 

Plenum  pressure  and  Plenum  pressure  error   (percent) 
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plotted  vs.  time  appear  in  Figs  5  and  6.  The  two  curves  are 
in  quite  good  general  agreement,  with  the  magnitude  of  error 
being  on  the  order  of  several  parts  per  thousand.  The 
linear  model  exhibits  a  greater  pressure  reduction  in  the 
initial  pressure  transient,  than  does  the  Loads  and  Motion 
results,  and  this  transient  is  more  heavily  damped  in  the 
linear  model.  One  explanation  for  this  difference  would 
appear  to  be  the  variable  leakage  area  of  the  Loads  and 
Motion  program.  Pitch  appears  to  couple  more  strongly  into 
pressure  in  the  linear  model;  however,  the  magnitude  of 
pitch  is  significantly  greater  in  the  linear  model  due  to 
the  absence  of  significant  damping  terms.  Finally,  the 
pressure  in  the  Loads  and  Motion  program  is  tending  toward  a 
slightly  higher  value,  further  supporting  the  supposition 
that  leakage  area  has  decreased  in  the  Loads  and  motion 
program. 

2 •   C-iG-i  Acceleration 

Figs  7  and  8  show  C.G.  acceleration  and  C.G. 
acceleration  error  vs.  time.  They  generally  exhibit  similar 
characteristics,  but  due  to  the  larger  initial  pressure 
reduction  of  the  linear  model,  there  is  a  more  pronounced 
overshoot  in  C.G.  acceleration  in  the  linear  model.  In 
addition,  the  stronger  pitch-pressure  coupling  in  the  linear 
model  also  translates  as  stronger  pitch-C.G.  acceleration 
coupling.  The  errors  appear  as  relatively  large 
percentages;  however,  these  are  the  differences  of  two  very 
small  numbers,  and  should  be  regarded  cautiously. 

3 .   Draft 

Draft   and  draft  error  vs.  time  appear  in  Figs  9  and 
10.   As  has  been   previously   verifed   by   the   Gerba/Ihaler 
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model,  dxaft  is  very  well  approximated  by  a  linear  model. 
For  this  particular  disturbance,  it  appears  to  be  slightly 
less  damped  and  tending  toward  a  lower  steady-state  value  in 
the  linear  mcdel.  Little  evidence  of  pitch  coupling  is 
apparent,  and  errors  are  on  the  order  of  parts  per  thousand. 

4.   Pitch  Rate 

Fitch  rate  and  pitch  rate  error  are  shown  plotted 
vs.  time  in  Figs  11  and  12.  The  linear  model  exhibits  a 
pitch  rate  cf  approximately  twice  the  magnitude  of  the  Loads 
and  Motion  program;  this  may  ostensibly  be  attributed  to  the 
sharper  pressure  transient  previously  observed;  with  the 
center  of  pressure  fprward  of  the  C.G.  in  both  models,  this 
would  induce  a  higher  initial  angular  acceleration  in  the 
linear  model.  In  addition,  the  linear  model  is 
significantly  less  damped  in  pitch;  therefore  this  larger 
magnitude  angular  acceleration  would  build  to  a  larger 
angular  velocity,  and  persist  for  a  longer  period  of  time  in 
the  linear  mcdel.  The  frequencies  compare  favorably,  5.04 
rad/sec  for  the  linear  model  and  5.28  for  the  Loads  and 
Mction  program. 

5«   Fitch  ftngle 

Figs  13  and  14  display  pitch  angle  in  degrees  and 
pitch  angle  error  vs.  time.  Both  show  basically  identical 
characteristics;  i.e.,  an  initial  pitch  down  induced  by  the 
pressure  drop,  and  an  exponential  decay  to  a  lower  pitch 
angle  induced  by  the  draft  change.  The  magnitude  of  pitch 
angle  change  reflects  the  larger  angular  velocity  of  the 
linear  mcdel,  and  the  difference  between  the  damping  of  the 
Loads  and  Motion  program  and  the  linear  model  is  quite 
dramatic.   The  magnitude  of  error,  however,  is  on  the   order 
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of    parts    per    hundred. 
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figure    5    -      PLENUM    PRESSURE    VS.    TIME 
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Figure  10  -   DRAFT  ERROR  VS.  TIME 
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Figure     11    -      PITCH    RATS, HAD/SEC    VS.    TIME 
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figure    12    -      PITCH    RATE    ERROR    VS.    TIME 
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Figure    13    -       PITCH    ANGLE,    DEGREES    VS.     TlflE 
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Figure  14  -   PITCH  ERROR  VS.  TIME 
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V.   FREQUENCY  DOMAIN  DEVELOPMENT 


This  chapter  develops  tha  techniques  for  frequency 
response  of  a  general  linear  system  by  complex  matrix 
inversion.  The  ocean  wave  equations  utilized  by  Booth  are 
reviewed  and  summarized,  and  the  forcing  functions  of  wave 
action  on  the  linearized  XR-3  model  are  developed  in  pitch 
and  heave.  Portions  of  Booth's  work  are  repeated  utilizing 
the  linearized  model,  and  the  results  contrasted  and 
compared . 


A.   LINEAR  SYSTEM  FREQUENCY  RESPONSE 


The  development  of  a  linear  system  permits  the 
utilization  of  various  techniques  for  the  development  of 
frequency  response  characteristics  which  are  vastly  superior 
in  CPU  time  in  simulation,  than  the  methods  to  which  one  is 
limited  in  non-linear  models.  The  definition  of  the  state 
matrix  of  the  previous  chapter  defines  the  state  transfer 
matrix: 


H(s) 


-1 
=|sI-A| *B (s) 


(7-1) 

and  if  the  system  is  forced  by  a  sinusoid,  s  tends  toward  jw 
as  the  transients  decay  and: 
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-1 

H  (jw)         =|  jm-A|*B(jw) 


(V-2) 


yielding  the  magnitude  and  phase  relationships; 
G(w)  =10*log*|H  (jw)/H(0)  | 


(V-3) 


9(w)  =tan|imag(H  ( jw)/real(H(  jw)  j 


(V-4) 


where  G(w)  is  the  gain,  in  decibels,  with  respect  to  some 
arbitrary  reference  H(0)  ,  and  9(«)  is  the  phase  angle  in 
the  complex  plane.  Thus  the  frequency  response  of  a  linear 
system  tc  a  given  sinusoidal  forcing  function  may  oe 
obtained  by  complex  matrix  inversion.  While  this  is  not  in 
itself  a  particularly  economical  use  of  CPU  time,  it  is 
still  vastly  superior  to  a  Fourier  transform.  This  is 
particularly  true  ).£  the  range  of  frequencies  is  wide  and 
variable-step  integrators  driven  to  small  step  sizes  which 
are  not  required  for  most  of  the  array  integration, 
finally,  it  is  not  necessary  to  allow  the  system  to  attain 
steady-state  AC  conditions  prior  to  computing  the  frequency 
response.  This  is  particularly  important  in  any  model  cf  the 
XR-3,  which  exhibits  very  light  pitch  damping. 


B.   WAVE  CHARACTERISTICS 

Booth  in  ref  £1J  provides  an  excellent  discussion  on  the 
characteristics  of  gravity-free  surface  waves  with  straight, 
infinitely  lcng,  parallel,  equally-spaced  wave   crests,   and 
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constant  wave  height.  Such  a  wave  advances  with  a  velocity 
C  which  is  often  referred  to  as  celerity,  in  order  to 
emphasize  the  fact  tnat  it  is  the  wave  front  which  advances, 
not  the  water  particles. 

The  regular  wave  may  thus  be  modeled: 

C  =G/v< 

l 

(V-5) 

Lambda        =2*Pi*C/w 

i 

2 
=2*Pi*G/w 

i 

(V-6) 

where   w  is   the   incident   wave   frequency,   Lambda   is  the 
i 

wavelength,  and  G  is  the   gravitational-  acceleration.    The 

wave  surface  is: 

A(x,t)        =A*cos (2*Pi*x/Lambda+w*t) 

0  i 

=  A*  (cos  (w*x/C) *sin  (w*t)  +sin  (w*x/C) *cos  (w*t) ) 
0      i  i        i  i 

(V-7) 

after  substitution  of  Equations  (V-5, 6) ,  and  the  appropriate 
trigonometric  identities.  The  frequency  experienced  by  the 
craft,  however,  is  the  encounter  frequency: 
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w  =  w*  (OV*cos  (Beta)  )  /C 

e  i 


2 
:w-w*V*cos  (Beta)/G 
i  i 


(7-8) 

where  V  is  the  craft  velocity  and  Beta  is  the  angle  of  craft 
heading  with  respect  to  the  advancing  wave  front,  Beta=0 
degrees  corresponding  to  the  craft  and  wave  front  travelling 
in  the  same  direction  (following  seas) . 


C.   WAVE  IOBCING  FUNCTIONS 


Wave   action   introduces  a  vertical  force  in  neave  and  a 
moment  about  the  center  of  gravity  in  pitch. 

It  is  assumed  that: 

1.  wave  generated  forces  and  moments  are  due  solely  tc  the 
net  additional  hull  volume  submerged  by  the  wave; 

2.  that  the  wave  does  not  exist  within  the  plenum. 

These  assumptions  appear  justified  due  to  the  small 
magnitude  of  waves  considered,  the  large  inertia  of  tne 
craft,  and  tte  filtering  action  of  the  bowseal. 

The  waveheight  as  a  function  of  the  x-location  along  the 
hull  and  time  t  may  be  described  by: 

Svht (x, t)      =Wvht*sin(w*t+2*Pi*x/Lambda) 

e 


=Wvht* [sin (w*t) *cos (x*w/C) 
L     e  i 
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♦  cos,(w*t)  *sin(X*v/cA 
e  i  J 


(  V-  9) 

after  substitution  from  equations  (V-5)  and  (V-6)  and  the 
appropriate  trigonometric  identities.  As  shown  by  Fig  15, 
the  differential  volume  of  both  sidewalls  immersed  by  the 
incident  wave  at  time  t  is: 

dV (x,t)       = (1/2)  *Wvht (x,t) *2*Ws*dx 

(  V-10) 

and  the  differential  heave  force  at  time  t  is: 
dHbuoyw(x,t)  =-Rho*G*dV  (x,  t) 

(  V-11) 

which  may,  after   appropriate   substitution   from   Equations 
(V-9,10) ,  be  integrated  analytically  from  stern  to  bow  as: 

Hfcuoyw  (t)     =-2*Rho*G*  (C/w)  *Wvht*sin(w*L/  (2*C) )  *sin  (w*t) 

i  i  e 

=HBUOYW*sin (w*t) 
e 

(  V-12) 


The  differential  moment  due  to  the   differential   volume 
submerged  by  the  incident  wave  at  time  t  is: 

dPbuoyw  (x,t)  =Rho*G*  (x-Xcg)  *dV  (x,t) 

=Rho*G*  (x-Xcg) *Wvht (x,t)  *Ws*dx 


(  V-13) 


which  integrates  from  stern  to  bow  as: 
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Pbuoyw  (t) 


•2*Rho*Ws*3*Wvht*cos  (w*t)  * 

e 


(  (C/w)  *sin  (w*L/(2*C)  ) 
i  i 


-  (C*L/(2*w)  ) *cos  (w*L/  (2*C) ) ) 

i  i 


+Xcg*Hbuoyw (t) 

=  PBUOYW*cos (w*t)  +Xcg*HBUOYW*sin (w*t) 
e  e 


(    V-1U) 

A  particular  note  of  caution  should  be   introduced   at   this 

point   with   regard   to   the   integrals   for  wetted  area  and 

wetted  moments   along   the   hull;   the   integrals   of   these 

functions   vanish   for  certain  w  .   Specifically,  for  heave, 

i 

Equation  (V-12)  vanishes  for  all  time  t  when: 


W*L/  (2*G) 

i 


=  n*Pi 


=  (2*n*Pi*G/L) 


1/2 


,11=1,2,3, 


(  V-15) 

The  pitch  integral  will  likewise  vanish  for  all  time  t  if, 
for  a  value  of  n  such  that  equation  (V-12)  is  zero,  the 
transcendental  equation  obtained  by  substitution  of  equation 
(V-15)  into  equation  (V-14) ,  and  equating  the  latter  to 
zero: 

2  1/2  2  1/2 

tan  (  (n*Pi*G*L/  (2*C))   )=  (2*n*Pi*G/  (L*C)  ) 

(  V-16) 
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is  also  satisfied.  Thus  the  frequency  response  of  the  craft 
should  exhibit  a  certain  ripple  effect  as  these  critical 
frequencies  are  passed. 

The  two  previously-defined   wave   integrals   HBUOYW   and 

PBQOYW   make  up  a  two-dimensional  B  vector  which  couples  the 

wave  force  onto  the  C.G.    acceleration  and  pitch  acceleration 

state   variables,  from  the  two  forcing  functions,  U=sin(w*t) 

1      e 

and  U=cos(w*t).   However,  it  is  convenient  to  redefine   this 
2     e 

relationship  in  terms  of  a  single  forcing  function  since  the 

sine  and  the  cosine  are  related  through  the  derivative.  In 
order  to  facilitate  this  relationship,  the  state  matrix  is 
augmented  with  a  dummy  state  variable: 

a  =  d/    (D"!=-w*  (PBUOYtf/Iyy)  *U 

10  /dt    2    e  1 

(  7-17) 

Thus  the  new  state  x   represents,   after   integration,   the 

10 

control   U  ,  multiplied  by  the  appropriate  B  (w)  vector  terms 
2         '  e 

of  PBUOYW/Iyy.   Note  that  if  a  time  domain  solution  of   wave 

forcing  action  were  desired,  a  new  initial  condition  would 
be  required,  representing  the  initial  wave  buoyant  moment; 
however,  that  does  not  enter  into  the  solution  in  the 
frequency  domain.  The  resulting  state  equations  are  shown 
in  Fig  16.  As  has  been  previously  discussed,  the  frequency 
response  may  now  be  computed  for  any  frequency: 


-1    I 

*B(w)  | 

e         (    e        e| 


G(jw)         =  |  jw*I-A|*B(w) i  | 


(  V-13) 
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where  G(jw)  is  the  column  vector  of   gains   for   each   state 
e 

variable. 


D.   FBEQOENCY  DOMAIN  SOLUTION 


The  previously-defined  state  equations   were   solved   in 

the   frequency    domain    utilizing   a   CSMP  III   computer 

simulation  which  appears  at   the   end   of   the   text.    This 

simulation   was   based   around  the  initialization  portion  of 

the  time  domain   solution   described   in   earlier   chapters; 

however,   rather   than   a  call  to  integration  in  the  Dynamic 

portion  cf  the  model,  the  CSMP  III   reserved   variable   TIME 

was   used   only   as   a   stepping   variable   to   generate   50 

frequencies   between   w  and   w,  the   specified   minimus   and 

0      max 

maximum  frequencies  respectively,  such  that: 

Delw*TIME 
w  =  w*e 

e  0 


(V-19) 

where   Delw  was  defined  such  that  when  TIME  =1.0,   w=w.   The 

e  max 

incident  wave  frequency  then  becomes,  for   ahead   seas,   the 

positive   solution   to  the   quadratic   Equation  (V-8)  ,  with 

cos  (Beta) =-1 .    The    influence    coefficients   about    the 

steady-state  operating  point  from  the  initialization  portion 

were  copied  cnto  a  complex  matrix  Ajw,  the   | jwI-A |   matrix, 

e 

and  the  E  vector  for  that  particular  incident  wave  frequency 

computed.  The  International  Mathematical  and  Statistical 
Litraries  (IMSL)  subroutine  LEQT2C,  complex  linear 
simultaneous  equation  solution,   was   used   to   compute   the 
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vector  of  gains: 

-1 
G(jw)  =|  jwI-A|*B  (w) 

e  e       e 


(7-20) 

and   the   resultant   absolute   values   of  G(jw)  were  plotted 

e 

against  the  encounter  frequency  w  as  decibels   with   resDect 

e 

to   the   unit   waveheight   for   heave,  heave  rate,  and  heave 

2 
acceleration,  and  with  respect  to  the   waveslope,   w/G,   for 

i 

pitch,   pitch   rate,  and  pitch  acceleration.   Note  that  this 

is  a  linear  system;  accordingly,  variations  in  waveheight  do 
not  affect  this  model's  response  characteristics,  as  was 
true  of  the  non- linear  Loads  and  Motions  program  in  ref  |1|. 


E.   COMPARISON  WITH  LOADS  AND  MOTIONS  PROGRAM 


Booth  in  ref  |1|  studied  the  frequency  response  of  the 
XR-3  Loads  and  Motion  program  in  response  to  seas  developed 
by  the  KAVZS  subroutine.  It  is  of  interest  to  compare  the 
linear  system  developed  in  this  study  with  the  results  of 
ref  |1|.  Duplicating  his  work,  three  runs  were  made 
representing  ahead  seas,  cross  seas  and  following  seas  at 
10,  20  and  30  kts  and  a  comparison  made  with  his  results. 
Merged  plots  of  pitch  and  heave  were  obtained  for  each  speed 
condition  for  encounter  frequencies  varying  from  1  to  25 
rad/sec.  In  addition,  Cal-Corap  plots  of  pitch  and  heave 
were  obtained  as  functions  of  the  ratio  of  craft  length  to 
wavelength.  Initial  conditions  and  craft  parameters 
selected  for  initialization  were  identical  to  those  selected 
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for   the   previous   chapter   with   regard    to   time-domain 
validation. 


1 .   Ahead  Seas  Run 

Fig  17  to  Fig  22  show  plots  of  linear  system 
response  in  heave  and  pitch  for  10,  20  and  30  knot  runs  into 
ahead  seas,  alternating  with  plots  reproduced  from  ref  £1J, 
representing  identical  runs  done  by  Booth  with  the  smallest 
wave  amplitudes.  These  were  selected  as  representing  the 
most  nearly-linear  mode  of  operation  of  the  Loads  and  Motion 
program  in  sea-state.  One  may  observe  the  general 
similarity  of  the  two  simulations,  with  the  speed-dependent 
resonant  peak  of  heave  in  the  linear  model  varying  from  an 
encounter  frequency  of  3.39  rad/sec  at  10  knots  ,  to  6.45 
rad/sec  at  30  kts,  while  the  non-linear  system  varies  from 
4.0  to  approximately  5.0  rad/sec.  The  frequency  interval 
of  the  non-linear  study  was  considerably  coarser  than  that 
utilized  in  the  linear  model,  1  rad/sec  vs.  less  than  .1 
rad/sec.  The  linear  system- is  significantly  less  responsive 
than  the  Loads  and  Motions  program,  by  11  dB  in  heave  and  6 
dB  in  pitch.  This  is  attributed  to  the  absence  of 
wave-plenum  chamber  interactions  with  the  plenum  pressure 
Pbtar,  and  indicates  that  the  presence  or  absence  of 
significant  wave  action  within  the  plenum  chamber  is  an  area 
of  considerable  impact  on  the  response  of  the  craft,  and 
should  be  included  in  future  studies. 

The  linear  system  is  quite  sharper  in  its  resonant 
peak  in  pitch,  as  should  be  expected  from  the  significantly 
reduced  pitch  damping  previously  demonstrated  by  the  linear 
model  in  the  time  domain  validation.  This  accounts  in  part 
for  the  tendency  of  the  linear  model  to  decrease  in  response 
slightly  from  an  encounter  frequency  of  1  rad/sec,  until 
just   before   the   onset   of  resonance,  unlike  the  Loads  and 
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Motions  data;  the  half-power  points  in  the  Loads  and  Motions 
frequency  resonance  are  more  widely  separated  than  are  those 
of  the  linear  system,  and  this  in  part  also  accounts  for 
this  difference. 

The  linear  system  shows  the  distinctly  periodic 
frequency  response  predicted,  with  a  sharp  null  occurring  at 

incident   wave   frequencies,   w,   of   3.14,   4.4,   and   5.38 

i 

rad/sec.    These   frequencies  correspond  to  a  value  of  n,  in 

eqn  (7-15),  cf  1,  2  and  3.  These  nulls  are  quite  sharp,  and 
are  not  as  well  defined  in  the  Loads  and  Motions  studies  due 
to  the  coarser  frequency  resolution,  but  are  nonetheless 
visible  in  all  runs. 

The  Loads  and  Motions  studies  also  demonstrate  an 
interesting  form  of  non-linear  response  with  regard  to 
resonant  pitch  frequency.  For  small  wave  amplitudes,  less 
than  .1  ft,  the  Loads  and  Motions  resonant  pitch  frequency 
varies  with  craft  speed  in  nearly  an  identical  manner  with 
the  linear  system,  with  the  slight  exception  of  the  30  kt 
study;  however,  the  only  available  data  for  that  run  was  at 
a  wave  amplitude  of  .2  ft.  For  wave  amplitudes  greater  than 
.1  ft,  the  resonant  peak  occurs  at  a  frequency  of  4  rad/sec, 
regardless  of  craft  speed.  A  tentative  explanation  for  this 
behavior  is  herein  advanced:  as  the  craft  resonant 
freguency  is  approached,  bowseal  gapping  commences  for  the 
larger  wave  amplitudes,  with  attendant  loss  of  plenum 
airmass;  the  loss  of  seal  force  and  plenum  pressure  cause 
the  craft  to  settle  back  into  the  water,  but  a  fixed  amount 
of  time  must  elapse  before  the  fans  can  replenish  the  plenum 
airmass  and  restore  normal  system  operation.  Thus  the  4 
rad/sec  period  observed  for  larger  wave  amplitudes 
represents  the  fastest  time  (.66  sec)  over  which  gapping 
can  occur  and  the  fans  sufficiently  replenish  the  airmass  to 


94 


perait  it  to  recur.  This  non-linear  response  is  somewhat 
analogous  tc  a  relaxation  oscillator,  and  represents  an 
interesting  area  for  future  study,  as  such  a  non-linearity 
might  easily  be  modeled  in  the  Nyquist  plane,  utilizing  the 
modified  Nyquist  stability  criteria  to  predict  this 
non-linear  behavior. 

a.   10  Kts  Ahead  Seas 

Figs  17  and  18  show  the  superimposed  heave  and 
pitch  response  of  the  linear  and  Loads  and  Motions  program, 
respectively,  for  an  ahead-seas  run  at  a  craft  velocity  of 
10  kts.  The  linear  jnodel  represents  the  response  of  the 
craft  in  heave  and  pitch  at  50  discrete  frequencies  ranging 
exponentially  from  1.0  to  25.0  rad/sec.  The  non-linear  data 
reproduced  from  ref  £ll  represents  8  discrete  frequencies 
over  that  same  range  at  1.0,  2.0,  3.0,  4.0,  5.0,  6.0,  8.0, 
and  10.0  rad/sec,  at  a  waveheight  of  .2  ft.  Although  this 
waveheight  was  the  smallest  used  by  Booth  for  runs  at  these 
speeds,  it  is  still  sufficiently  large  to  excite  the 
aforementioned  non-linear  response.  This  behavior  makes 
itself  manifest  as  the  4.0  rad/sec  resonance  in  pitch, 
compared  fcith  3.39  rad/sec  for  the  linear  model.  The  Loads 
and  Motions  program  shows,  for  small  weight  disturbances  at 
these  speeds,  a  natural  pitch  period  of  1.64  sees,  or  3.7 
rad/sec.  As  previously  noted,  the  magnitude  of  the  linear 
system  pitch  response  agrees  fairly  well,  within  3  dE,  of 
the  non-linear  Loads  and  Motions  data,  with  a  sharper  peak 
due  to  the  lighter  linear  damping;  the  heave  magnitudes 
differ  considerably,  by  -15  dB,  due  to  the  absence  of  plenum 
waves  in  the  linear  system. 

kith  the  exception  of  the  heave  magnitude 
difference  and  the  shift  in  resonance  due  to  the  ncn-linear 
behavior,   the  shape  of  the  heave  response  curves  of  the  two 
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models  compare  quite  favorably.  A  shallow  null,  of 
undetermined  origin;,  appears  in  both  heave  curves  at  3.0 
rad/sec,  just  slightly  before  the  onset  of  resonance.  A 
very  sharp  null  occurs  at  8.36  rad/sec  in  bcth  models;  this 
is  at  an  incident  wave  frequency  of  3.15  rad/sec,  or  a 
wavelength  of  20  ft.  As  predicted  by  Equations  (¥-15,16) , 
further  nulls  occur  at  a  wave  incident  frequency  of  4.46  and 
5.58  rad/sec,  corresponding  to  wavelengths  of  10  and  6.6 
feet  respectively,  but  the  Loads  and  Motion  encounter 
frequencies  did  not  range  sufficiently  high  to  observe  them. 

Pitch  also  shows  general  agreement;  but  Booth's 
Loads  and  Motion  data  shows  a  null  at  6.0  rad/sec,  vs.  11.54 
rad/sec  for  the  linear  system.  This  null  may,  however,  be 
related  to  the  non-linear  response  previously  described, 
representing  a  freguency  at  which  the  phase  relationship 
between  plenum  airmass  replenishment  and  pitch  buoyancy  due 
to  waves  is  such  that  cancellation  occurs.  It  is 
suspiciously  located  at  a  3/2  multiple  of  the  non-linear 
response  frequency  of  4.0  rad/sec.  In  addition,  Booth  noted 
this  identical  discrepancy  in  comparison  with  the  linearized 
Aerojet  XE-3  model;  the  Aerojet  model  also  predicted  a  null 
in  the  vicinity  of  12  rad/sec,  which  Booth  was  unable  to 
observe  with  the  Loads  and  Motions  program,  but  which  is 
distinctly  visible  in  the  linear  model  at  11.54  rad/sec. 

b.   20  Kts  Ahead  Seas 

Figs  19  and  20  show  the  linear  and  Loads  and 
Motions  heave  and  pitch  frequency  response  for  a  20  let  run 
into  ahead  seas.  The  Loads  and  Motions  data  is  for  a  wave 
height  of  .1  ft  and  %s  more  nearly  linear  in  behavior  than 
previously  observed  at  10  kts.  The  range  of  frequencies 
spanned  by  the  linear  system  is  the  same  as  before;  tne 
Loads   and   Motion  data  includes  two  additional  points  at  15 
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and  20  rad/sec  in  heave  response.  The  linear  pitch  response 
is  smaller  than  that  observed  at  10  kts,  being  about  6  dB 
les.s  than  the  non-linear  data,  while  the  heave  response 
remains  consistently  less  than  the  non-linear  data  by  12  d3 
as  at  10  kts. 

The  shape  of  the  two  sets  of  frequency  response 
data  is  in  guite  good  agreement  for  these  conditions,  as  the 
Loads  and  Motions  is  apparently  in  a  nearly-linear  mode  of 
operation.  The  pitch  data  shows  a  distinct  flattening  at 
the  peak  of  resonance  in  both  models,  although  this  occurs 
somewhat  earlier  in  the  non-linear  data  (4.0  rad/sec)  than 
in  the  linear  data (4.68  rad/sec),  due  to  the  more 
sharply-peaked  pitch  response  of  the  latter.  Both  models 
resonate  at  5.0  rad/sec,  and  fall  off  above  resonance  with 
identical  rcll-offs  of  12  dB  per  octave,  after  the  initially 
steeper  roll-off  of  the  linear  system  from  peak  resonance. 
Heave  shows  the  characteristic  null  in  both  models  at  an 
encounter  freguency  of  14  rad/sec;  this  is  a  wave  freguency 
of  3. 15  rad/sec,  or  a  wavelength  of  20  ft.  The  linear 
system  predicts  a  pitch  null  at  20.0  rad/sec,  but  the  Loads 
and  Motions  data  does  not  extend  sufficiently  high  to 
observe  this. 

c.   30  Kts  Ahead  Seas 

figs  21  and  22  show  the  linear  and  ncn- linear 
freguency  response  for  a  30  kt  ahead-seas  run.  As  in  the  10 
kt  run,  the  smallest  wave  amplitude  utilized  in  Booth's 
study  was  .2  ft;  however,  it  is  considered  that  tne 
nose-down  planing  moment  which  is  quite  strong  at  high 
speeds  restrains  bowseal  gapping,  and  the  Loads  and  Motions 
data  represents  conditions  for  the  non-linear  response  which 
is  just  en  the  verge  of  becoming  established.  This  is 
supported   by  the   observation   that  pitch  resonance  occurs 
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between  the  frequency  previously  observed  as  the 
characteristic  non-linear  response,  and  linear  resonance, 
between  4.0  and  5.0  rad/sec;  furthermore,  the  heave 
resonance  is  weakly,  peaked  at  a  slightly  higher  frequency 
than  pitch  response,  approximately  5.0  rad/sec,  still  closer 
to  the  predicted  linear  resonance  at  6.46  rad/sec. 

The  linear  model  shows  still  lower  magnitudes  of 
pitch  response,  6  to  8  dB  below  the  Loads  and  Motion  data. 
The  heave  magnitudes  remain  consistently  12  d3  below  the 
Loads  and  Motion  data. 

Fitch  resonates  at  6.46  rad/sec  in  the  linear 
system,  and  between  4.0  and  5.0  rad/sec  in  the  non- linear 
Loads  and  Motions  data.  Roll-off  in  pitch  remains  the  same 
above  resonance,  approximately  12  dB  per  octave  in  both 
models,  but  heave  response  in  both  models  falls  cff  more 
shallowly,  approximately  10  dB  per  octave  in  both  models. 
The  heave  null  due  tp  the  20  ft  wavelength  is  still  apparent 
in  the  linear  model,  representing  an  encounter  frequency  of 
20.0  rad/sec,  but  the  Loads  and  Motions  data  did  not  extend 
sufficiently  high  in  frequency  for  this  null  to  be  observed. 

d.   Frequency  Response  vs.  Wavelength 

Pigs  23  and  24  are  plots  of  heave  response,  vs. 
the  ratio  of  craft  length  to  wavelength  for  the  linear 
system,  and  vs.  the  ratio  of  wavelength  to  craft  length  for 
the  Loads  and  Motion  data  reproduced  from  ref  £  1 J  -  Both 
plots  are  superimposed  with  the  data  from  all  three  velocity 
runs  of  10,  20,  and  30  kts.  The  difference  in  scale 
selected  for  plotting  in  the  linearized  model  precludes 
crowding  at  the  low  end  of  the  x-axis,  since  the  bulk  cf  the 
frequencies  studied,  3.15  rad/sec  and  above,  are  shorter 
than  craft  length. 
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The  linear  system  heave  response  shows  distinct 
nulls  for  waves  which  are  integer  multiples  or  submultiples 
of  the  craft  length,  and  may  be  observed  at  L/Lambda= 
1/4,1-0,2.0,  and  3.0;  less  distinct  nodes  may  be  observed  at 
L/lambda=  3/2,  5/2.r  and  7/2.  The  corresponding  nulls  at 
Lambda/L=4 . 0#1.0,  and  .5  are  visible  in  Fig  24,  but  scale 
crowding  and  the  coarser  frequency  resolution  utilized  by 
Booth  render  these  less  obvious.  The  sharp  node  at 
L/Lambda=.3  corresponding  to  a  wave  frequency  of  1.75 
rad/secf  is  net  the  result  of  the  periodic  hull  integrals; 
this  wave  frequency  corresponds  to  the  encounter  frequency 
of  3.4  to  6.C  rad/sec,  and  is  the  craft's  natural  pitch 
frequency.  It  may  te  observed,  also  as  a  node,  in  the  pitch 
response.  The  natural  pitch  frequency  variation  with  speed 
is  obscured  in  this  plot,  since  as  the  speed  increases,  the 
natural  pitch  frequency  increases;  but  this  is  an  encounter 
frequency,  and  the  wave  incident  frequency  which  produces 
this  encounter  frequency  decreases  with  increasing  speed. 
Therefore,  the  two  effects  nearly  cancel. 

Figs  25  and  26  show  pitch  response  for  the 
linear  system  and  the  non-linear  loads  and  Motion  data  of 
ref  £  1J  ,  plotted  as  was  heave.  Note  that  the  nulls  and 
nodes  are  roughly  180°  out  of  phase,  a  node  in  pitch 
corresponding  to  null  in  heave  and  vice-versa. 

2-   Cross-seas  Runs 


Figs  27  to  31  are  presentations  of  heave  and  pitch 
frequency  response  of  the  linear  system  and  the  Loads  and 
Motions  data  from  ref  £l J  at  20  and  30  kts  under   cross-seas 

conditions    (w=w) .     Comments    made    regarding   the   the 

e  i 

ahead-seas  runs   are   generally   applicable   here,   but   the 
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comparison  of  the  two  systems  is  in  no  manner  nearly  as 
striking  as  thar  for  the  ahead-seas  runs.  This  is 
attributed  to  the  significant  effects  of  roll  and  yaw  modes 
on  the  craft  dynamics  under  these  conditions,  which  are 
completely  neglected  in  this  linear  model. 

3 •   Following  Seas  Huns 

Booth  acknowledged  little  agreement  between  his  work 
with  the  Loads  and  Motion  program  and  the  Aerojet  XR-3 
linear  model  freguency  response  curves;  this  is  apparently 
due  to  an  error  in  calculating  the  encounter  freguency. 
Higher  freguency  waves  propagate  at  a  slower  velocity,  so 
that  for  a  craft  traveling   at   a   given   velocity,   as   the 

incident   wave   freguency   w  increases,   a  freguency  will  be 

i 

attained  above  which  the   craft   will   actually   outrun   the 

following  wave,  and  plow  into  such  a  wave  bows-on.  Thus  the 
craft  will  actually  experience  an  ahead-sea  condition   above 

that   wave   freguency.    The   w  corresponding  to  the  maximum 

i 

positive  w  that  the  craft  may  experience  for  following   seas 
(cos  (Beta)  =  1 .0)  is  given  by: 


2 

d/dw{w-V*w/G)  =0.0 
l  i    i 


(V-21) 
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w  =  G/(2*V) 

i 


(V-22) 


Substituting    Equation     (V-22)     into     (V-8)     yields: 

2 
v  =  G/  (2*V)-  (V/G)  *  (G/  (2*7)  ) 

emax 

=G/  (4*V) 


(V-23) 

Incident  wave  frequencies  above  or  below  the  w  corresponding 

i 

to  w      result   in   a   lower   w,   and   if   w  is   increased 
emax  e  i 

sufficiently,  past  w=G/V,  w  will  become  negative,  indicating 

i      e 

that  the  craft  is  overrunning  the  wave.  Accordingly,  a  great 

deal  of  judiciousness  must  be  utilized  when  interpreting  the 

results  of   a   following   seas   run,   when   plotted   against 

encounter  frequency  w. 

e 


Accordingly,  the  following-seas   frequency   spectrum 

was   generated   in   a  manner  different  from  that  utilized  on 

previous  runs.   W  was  varied  exponentially  over  the  range  of 

i 

1   to   25   rad/sec,   and  the  corresponding  w  was  calculated. 

Hence,  the  results  are   not   comparable   directly   with   the 

results   of  Booth,  but  are  presented  in  Fig  32  to  Fig  33  for 
information  only,  at  v=20  lets. 
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Figure    15    -       HAVE    DIFFERENTIAL     fORCES    AND    MOMENTS 
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Figure    17    -       LINEAR    SYSTEM    FREQUENCY    RESPONSE, 
10    KTS,    AHEAD    SEAS 
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V=10  lets,    A=0.2  ft, 
Ahead  Seas 
Q  =  Heave  Frequency  Response 

O  =  Pitch  Frequency  Response 

Figure  18  -   LOADS  AND  MOTION  FE2QUENCY  RESPONSE, 

10  KTS,  AHEAD  SEAS 
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Figure    19    -       LINEAR    SYSTEM    FREQUENCY    RESPONSE, 
20    KTS,    AHEAD    SEAS 
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Figure    20    -       LOADS    AND    MOTION    FBEQUENCY    RESPONSE, 

20    KTS,     AHEAD    SEAS 
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figure    21    -       LINEAR    SYSTEM    FREQUENCY    RESPONSE, 
30    KTS,    AHEAD    SEAS 
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Figure  22  - 


V=30  kts,   A=0.2  ft, 
Ahead  Seas 

Heave  Frequency  Response 
Pitch  Frequency  Response 

LOADS    AND    MOTION    FREQUENCY    RESPONSE, 
30    KTS,     AHEAD    SEAS 
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Figure  23  -   LINEAR  HEAVE  FREQUENCY  RESPONSE  VS.  L/LAMBDA 
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Figure  24  -   LOADS  AND  MOTION  HEAVE  RESPONSE  VS.     LAHBDA/L 
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Figure  25  -   LINEAR  PITCH  FREQUENCY  RESPONSE  VS.  L/LAM3DA 
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Figure  26  -   LOADS  AND  MOTION  PITCH  SESPONSE  VS.  LAKBDA/L 
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Figure    27    - 


LINEAR    SYSTEM    FREQUENCY    RESPONSE, 
20    KTS,    CROSS    SEAS 
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Figure  23  -   LOADS  AND  MOTIONS  HEAVE  RESPONSE, 
20  AND  30  KTS,  CROSS  SEAS 
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figure  29  - 


V=20  ktS,   A=0.2  ft, 

Abeam  Seas 
=  Pitch  Frequency  Response 
=  Roll  Frequency  Response 

LOADS  AND  MOTION  PITCH  RESPONSE, 
20  KTS,  CROSS  SEAS 
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Figure    30    - 


LINEAR    SYSTEM    FREQUENCY    RESPONSE, 
30    KTS,    CROSS    SEAS 
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Figure  31  - 


V=30  lets,   A=0.2  ft, 

Abeam  Seas 

Pitch  Frequency  Response 

Roll  Frequency  Response 

LOADS  AND  MOTION  PITCH  RESPONSE, 

30  KTS,  CROSS  SEAS 
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Figure    32    -       LINEAR    SYSTEM    FREQUENCY    RESPONSE, 
20    KTS,     FOLLOWING' SEAS 
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VI.   FREQUENCY  DOMAIN  REDUCTION  QI    XR-3  DATA 


This  chapter  develops  the  techniques  for  fast  efficient 
reduction  of  XR-3  testcraft  data  from  the  time  domain  to  the 
frequency  domain.  Some  more  prominent  noise  sources  are 
identified,  and  the  mode  of  noise  transmission  is 
hypothesized.  A  very  limited  amount  of  pitch  data  is 
collected  and  closely  examined  using  these  techniques  and 
compared  with  the  theoretical  development  of  the  preceding 
chapters. 


DATA  ACQUISITION  SYSTEM 


The  data  acquisixicn  system  consists  of  the  standard 
semi-per manently  installed  Pemco  tape  recorder.  This  is  a 
14-channel  FM  tape  system,  with  output  0-1  volt  positive.  A 
19-channel  patch  panel  in  the  tape  recorder  compartment, 
located  just  aft  of  the  cockpit,  permits  selection  of  the 
data  to  be  recorded; voice  edge  track  recording  directly  from 
the  craft's  radio/ICS  system  permits  relatively  simple 
documentation  of  experimental  runs  to  be  maintained. 


B.   FREQUENCY  DOMAIN  DATA  REDUCTION  SYSTEM 


The  £io- Engineering  Laboratory,  Electrical  Engineering 
Department,  NPS,  Monterey,  CA,  has  a  very  powerful  system 
for  spectral  analysis  of  analog  data.  This  system, 
consisting    cf    a    PDP-11    Mini-computer   and   Time-Data 
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Corporation  Superpanel-50  Sp.ectrum  Analyzer,  was  utilized 
for  time-to-frequency  domain  data  reduction  of  the  analog 
tape  data.  This  system  is  capable  of  performing  a  Fast 
Fourier  Transform  (FFT)  on  analog  data  of  magnitudes  of  +5V, 
from  a  maximum  upper  frequency  of  10  Hz  to  10  MHz,  for  up  to 
8094  discrete  time  samples.  However,  caution  must  be 
exercised  tc  avoid  memory  overflow  and/or  program  interrupt 
of  the  PDP-11  for  samples  in  excess  of  2048.  In  addition, 
this  entire  process  may  be  repeated  up  to  2048  times,  and 
the  resultant  Fourier  coefficients  averaged,  for  further 
reduction  of  transient  frequency  components.  The  resultant 
transformed  data  may  be  output  as  Fourier  coefficients, 
power  spectral  coefficients,  and  as  an  amplitude  histogram; 
output  devices  include  a  storage  CRT  and  a  Hewlett-Packard 
70C4-B  X-Y  plotter  which  will  reproduce  the  display  on  the 
CRT  on  ccmraand  from  the  panel.  In  addition  the  time  domain 
data  being  analyzed  may  be  observed  on  a  time-base  storage 
oscilloscope,  and  monitored  throughout  the  transformation 
process.  Processing  times  vary,  of  course,  with  upper 
frequency,  cumber  of  samples,  and  number  of  averaging 
"windows",  fiom  less  than  a  second  to  several  minutes.  For 
further  information  see  ref  [8 J ,  the  operating  instructions 
for  the  Superpanel-50 .  Fig  34  shows  the  details  of  the 
setup  for  data  reduction  to  the  frequency  domain. 


C.   FREQUENCY  DOMAIN  ANALYSIS 


XR-3  data  was  collected  on  an  experimental  run  conducted 
20  January,  1977  at  Lake  San  Antonio  for  the  purpose  of 
performing  frequency  domain  analysis,  using  the 
previously-described  FFT  techniques.  Two  analyses  were 
performed : 

1.   a   broad-band   noise   analysis,  examining  the  frequency 
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spectra  of  all  channels  up  to  200  Hz,  as  various 
suspected  noise  sources  were  brought  up  and  shut  down, 
or  varied  in  power  setting;  the  purpose  of  this 
experiment  was  to  identify  offending  noise  sources  by 
their  characteristic  signature  in  the  frequency  domain; 

2.  a  lcw-f reguency  analysis  of  craft  characteristics  in 
pitch  and  pitch  rate,  for  15  kt  runs  into,  across,  and 
with  the  prevailing  wave  directions.  The  frequency 
spectra  of  these  channels  were  analyzed  from  .025  to 
25  Hz,  for  the  purpose  cf  comparison  with  ref  £lj  and 
the  work  cf  the  preceding  chapters. 

1  •   Wideband  Noi-sa  Analysis 

Noise  data  from  the  XR-3  was  analyzed  over  a 
freguency  range  up  to  200  Hz  and  512  samples,  averaged  over 
8  windows;  this  gave  a  total  sample  time  of  approximately  30 
seconds.  A  Butterworth  filter  was  used  to  reject  signals 
above  200  Hz  to  prevent  aliasing  of  higher  freguency 
components  into  the  lower  range  of  interest  due  to 
discretization.  In  addition,  a  battery/rh-eostat  combination 
was  used  to  bias  out  the  1/2  volt  offset  of  the  F3  tape, 
since  ,  as  previously  noted,  the  FM  tape  output  range  is  0-1 
vclt  positive,  while  the  input  to  the  Analog-to-Digital 
converter  of  the  Super panel-50  is  +5V;  this  considerably 
improves  amplitude  resolution  of  the  transformed  data. 

A  representative  output  of  freguency  spectra  appears 
in  Fig  35;  the  remainder  may  be  examined  in  Appendix  A.  The 
data  from  all  channels  are  basically  identical  in  nature, 
and  too  repetitive  to  include  in  the  main  text.  As  nay  be 
expected,  the  craft  main  engines  and  the  60-Hz  power  supply 
contribute  tte  major  share  of  noise.  As  power  is  varied,  the 
main  engine  noise  may  be  observed  to  vary  from  20  Hz  tc   100 
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Hz,  while  the  power  supply  contributes  significant  noise  at 
55-60  Hz  and  its  harmonics.  Following  these  sources  in 
order  of  magnitude  are  the  lift  fans,  which  may  be  observed 
in  the  vicinity  of  40  Hz. 

An  attempt  was  made  to  ascertain  the  mode  of 
transmission  of  this  noise,  specifically,  whether  by 
electrical  means,  or  by  mechanical  vibration  of  the  sensor 
package.  Accordingly,  the  gyro  package  containing  pitch, 
roll,  and  yaw,  and  pitch  rate,  roll  rate,  and  yaw  rate  gyros 
was  removed  and  placed  on  a  Ling  shaker  table,  and  exposed 
to  vibrations  of  1/2  to  1  G,  over  a  frequency  range  of  5-200 
Hz.  The  RMS  shaker  table  displacement  was  measured  with  a 
calibrated  Bentley  accelerometer ,  and  compared  with  the  RMS 
voltage  output  of  the  individual  gyros.  It  was  hoped  to 
obtain  frequency  response  curves  of  the  gyro  packages,  but 
they  seemed  quite  insensitive  to  vibrations  in  this 
frequency  range,  and  the  frequency  curves  showed  significant 
scatter  due  to  the  low  signal-to-noise  ratio.  They  were  not 
considered  worthy  of  inclusion. 

It  was  concluded,  on  the  basis  of  this  experiment, 
that  mechanical  vibration  could  not  transmit  such 
overwhelming  amounts  of  noise  as  was  observed  being  injected 
into  the  XR-3  data  collection  system.  The  medium  of 
transmission  is  thus  electrical;  this  conclusion  is  further 
supported  by  the  observation  that  the  eight  engines 
installed  aboard  the  XR-3  have  no  ignition  noise 
suppressors;  likewise,  the  60  Hz  power  supply  is  unfiltered. 
Furthermore,  while  bonding  straps  run  the  length  of  the 
craft,  the  possibility  of  extensive  AC  ground  loops  is 
always  present  in  a  non-conducting  hull,  such  as  the  XR-3. 
These  areas  should  receive  some  close  scrutiny  in  the 
future,  for  a  tenfold  reduction  in  noise  injection  into  the 
data  collection  system  should  be  possible,  by  means  of 
ignition  noise  suppression,  and  the  installation   of   .1   uF 
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capacitors   at  each  AC  outlet,  at  a  relatively  insignificant 
cost  in  time  and  money. 

2 •   Low,-Ereguency  Craft  Dynamics 

Three  runs  were  made  during  the  aforementioned  data 
collection  runs,  in  which  pitch  and  pitch  rate  were  recorded 
for  15  kt  runs  into,  across,  and  with  the  prevailing  seas. 
Qualitative  sea  conditions  at  the  time  of  observation  were 
noted  as  wave  heights  of  3  to  6  inches,  with  a  wavelength  of 
approximately  6  feet_;  this  wavelength  represents  an  incident 
wave  freguency  of  approximately  2.2  rad/sec,  corresponding 
to  encounter  freguencies  of  6.2,  2.2,  and  -1.8  rad/sec,  for 
ahead,  cress,  and  following  seas  at  15  kts.  Data  was 
reduced  in  a  manner  similar  to  that  used  for  wideband  noise, 
with  the  exception  that  a  maximum  freguency  of  25  Hz,  1024 
samples  with  no  averaging  was  utilized.  »  A  difficulty 
developed  in  that  a  Butterworth  filter  with  cut-off  below 
100  Hz  was  not  readily  available  at  the  time.  However,  when 
cne  is  constrained  from  lowering  the  bridge,  it  may 
nevertheless  be  possible  to  raise  the  river.  In  this 
particular  case,  the  tape  was  run  at  four  times  recorded 
speed,  or  7  1/2  inches  per  second.  The  FFT  was  set  for  100 
Hz,  as  was  the  Butterworth  filter,  and  the  desired 
resolution  of  102Q  samples  selected.  This  procedure  had  the 
beneficial  side-effect  of  reducing  the  real  elapsed 
processing  time  of  the  Super panel-50/PDP- 1  1  combination  by  a 
factor  of  four.  The  resultant  transformed  data  was 
scale-expanded  after  computation,  so  that  only  the  portion 
from  0-5  Hz  (0-1.25  actual  Hz)  was  plotted. 

These  admittedly  unorthodox  and  awkward  procedures 
were  necessitated,  not  only  by  the  lack  of  suitable 
filtering,  but  also  by  difficulties  experienced  in  obtaining 
proper  operation  of  the  PDP-11  for  very  low   freguencies   at 
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that  particular  time.  Long  sample  times  combined  with 
frequent  inexplicable  program  interrupts  rendered  this 
particular  session  most  frustrating.  However,  ventilation 
and  overheating  problems  were  apparently  the  major  source  of 
difficulties  experienced  with  this  particular  exercise,  and 
future  data  reductions  in  the  frequency  domain  should 
reguire  ncne  of  these  awkward  circumlocutions. 

a.   Fitch  and  Pitch  Bate-Ahead  Seas  Bun 

figs  36  and  37,  the  reduced  FFT  of  pitch  and 
pitch  rate  for  the  15  Jet  ahead-seas  run,  may  be  compared 
with  Figs  38  and  39-,  the  linear  system  theoretical  response 
in  pitch  and  pitch  rate.  The  linear  system  was  plott€d  as 
dB  with  respect  to  waveslope  logarithmically  against 
encounter  frequency,  on  a  range  of  .1  to  7.85  rad/sec.  The 
observ-ed  data  are  relative  magnitudes  plotted  linearly  vs. 
encounter  frequency,  on  a  range  of  0  to  7.85  rad/sec.  The 
unusual  scale  was  necessitated  for  the  observed  data,  to 
obtain  the  conversion  to  rad/sec;  in  Hz,  the  scale  is  .125 
Hz/in.  A  qualitative  comparison  of  the  two  plots  shows  a 
rather  remarkable  similarity;  the  FFT  of  pitch  shews  a 
falling  off  in  response  from  the  very  low  frequency  end  to  a 
local  minima  at  2.0  rad/sec  with  the  onset  of  resonance 
occurring  in  the  vicinity  of  2.75  rad/sec;  resonance  peaks 
at  3.9  rad/sec  and  is  completed  at  approximately  5.0 
rad/sec.  Two  small  nodes  may  be  observed  at  2.36  and  6.23 
rad/sec,  the  latter  corresponding  to  the  encounter  frequency 
of  the  qualitatively-observed  wavelength  at  the  time  of  this 
run;  this  was,  as  previously  noted,  an  incident  wave 
frequency  of  2.2  rad/sec. 

The  linear  model  replicates  the  observed  data  in 
pitch  almost  precisely,  with  the  exception  of  the  two  nodes 
at  2.2  and  6.28  rad/sec,  due   to   the   dominant   local   wave 
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amplitudes  in  the  observed  data. 

As  may  be  expected,  pitch  rate  data  is  somewhat 
more  noisy  than  pitch  data,  although  not  excessively  so. 
ctserved  pitch  rate  shows  an  unusual  double  peak  in 
resonance,  which  also  occurs  in  the  other  two 
sea-conditions.  The  first  peak  occurs  at  4.35  rad/sec, 
slightly  higher  in  frequency  than  the  observed  pitch  data, 
while  the  second  peak  occurs  at  5.25  rad/sec  and  is  of 
unexplained  origin.  This  peak  may  be  indicative  of  the 
onset  of  the  non-linear  behavior  discussed  in  the  preceding 
chapter,  but  that  is  only  a  speculative  observation. 

The  linear  system  pitch  rate  data  shews  a 
sharper  resonance  than  did  linear  pitch,  the  former  being 
markedly  peaked  at  4.26  rad/sec,  while  the  latter  was 
relatively  flat  from  3.9  to  4.26  rad/sec.  The  agreement  in 
the  sharpness  of  pitch  rate  response  between  the  observed 
data  and  the  linear  system  is  quite  good.  Resonant  rise 
also  occurs  earlier  in  the  linear  model,  at  approximately 
2.0  rad/sec;  but  there  is  no  indication  of  the  peak  at  5.25 
rad/sec  which  occurs  in  the  observed  data. 

No   significant   nodes   or   nulls   occur  in  this 

frequency  range;  the  pitch  node  at  n=1 (  w=3.15  rad/sec)  does 

i 

not  occur  until  an  encounter  frequency  of  9.75  rad/sec. 

Overall,  the  observed  and  linear  data  compare 
remarkably  well  under  these  conditions. 


b.   Pitch  and  Pitch  Rate-Cross  Seas  Run 

Despite   the   fact   that   cross-seas  runs  induce 
modes  of  motion  in  roll  and  yaw  that  are  unaccounted  for   in 
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this   linear  model,  it  was  considered  of  interest  to  compare 

observed  cross-seas  data,  Figs  40  and  41,  with  Figs   42   and 

43,  the  linear  data  for  w  =  w. 

e  i 

The  observed  pitch  FFT,  Fig  40,  is  more 
distinctly  resonant  at  4.0  rad/sec  under  these  sea 
conditions  than  in  the  ahead-sea  run;response  decreases 
continuously  from  the  very  low  frequency  end  until  almost 
the  onset  cf  resonance  at  approximately  3.5  rad/sec. 
Resonance  occurs  at  4.0  rad/sec,  and  a  distinct  null  occurs 
at  5.25  rad/sec.  The  observed  predominant  wavelength,  which 
as  noted  previously  corresponds  to  an  incident  wave 
freguency  cf  2.2  rad/sec,  is  distinctly  visible  at  that 
freguency. 

Comparison  with  the  linear  system  pitch  response 
curves  cf  Fig  42  is  as  striking  as  the  ahead-seas  runs,  the 
more  so  because  of  the  modes  of  motion  induced  in  the  craft 
which  are  not  in  the  linear  model.  As  in  the  FFT  data, 
linear  pitch  response  decreases  until  almost  the  onset  of 
resonance,  with  a  pronounced  null  at  3.5  rad/sec.  This  null 
does  not  appear  in  the  observed  pitch  data,  but  is 
distinctly  visible  in  the  observed  pitch  rate  data.  linear 
system  resonance  is  sharper  under  these  sea  conditions,  as 
is  the  observed  data,  and  occurs  at  4.26  rad/sec.  Two 
high-frequency  nulls  at  5.25  and  7.0  rad/sec,  and  a  node  at 
6.0  rad/sec  ,  are  distinctly  visible  in  the  observed  data, 
although  less  prominently. 

The  observed  pitch  rate  data  shows  a  sharper 
rcll-off  from  the  low  freguency  end  of  the  spectrum,  and  a 
flatter  response  thereafter  until  the  onset  of  resonance, 
than  does  pitch.  The  distinct  double  peak  previously  noted 
in  the  ahead-sea  run,  is  distinctly  visible  in  this  run 
also. 
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Ccmparispn  of  the  pitch  rate  FFT,  Fig  41 ,  with 
the  linear  system  pitch  rate  response  of  Fig  43  shows  the 
same  good  agreement  as  did  the  pitch  data;  the  low  freguency 
roll-off  characteristics  are  very  much  the  same,  and  the 
null  predicted  at  3.5  rad/sec  by  the  linear  model  is 
prominently  visible  in  the  FFT.  However,  the  higher 
freguency  nodes  and  nulls  are  much  less  distinct. 

c.   Pitch  and  Pitch  Rate-Following  Seas 

Due  to  the  difficulties  previously  experienced 
in  interpreting  results  of  following  seas  runs  with  Booth's 
work,  only  a  most  general  comparison  of  observed  pitch  and 
pitch  rate  FFT's,  shown  in  Figs  44  and  45,  and  the  linear 
system  responses,  shewn  in  Figs  46  and  47,  will  be  made. 
Note  that  due  to  the  technique  utilized  in  generating 
encounter  frequencies  in  the  linear  model,  the  range  of 
frequencies  plotted  in  the  linear  model  is  much  greater,  (.3 
rad/sec  to  -40  rad/sec) ,  resulting  in  a  much  coarser 
frequency  resolution  in  the  range  of  interest.  Negative 
encounter  frequencies,  as  noted  in  the  preceding  chapter, 
represent  a  modified  ahead-sea  condition. 

The  linear  system  predicts  a  sharp  roll-off  in 
pitch  response  from  the  very  low  frequency  end  to  about  -.25 
rad/sec,  and  a  relatively  flat  response  from  there  until  the 
onset  of  resonance,  at  about  -2.0  rad/sec.  Resonance  occurs 
quite  sharply  at  -4.11  rad/sec.  No  nodes  or  nulls  in  this 
frequency  range  are  predicted  by  the  linear  model.  The 
observed  pitch  data  is  in  general  agreement  with  this;  in 
addition  the  predominant  wavelength,  which  occurs  under 
these  sea  conditions  at  an  encounter  frequency  of  -1.6 
rad/sec,  may  be  observed  in  the  FFT. 

The  linear  system  and  observed  data  do  not  agree 
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well  in  pitch  rate,  the  linear  system  predicting  an 
increasing  response  from  the  low-frequency  end  until 
resonance,  while  observed  data  shows,  as  should  be  expected, 
a  low-frequency  roll-off  to  a  relatively  flat  mid-range 
response.  No  plausible  explanation  may  be  advanced  for  this 
variance  at  this  time;  however,  considerable  difficulty  has 
been  observed  by  both  Booth  and  this  author  in  correctly 
interpreting  following  seas  results. 

The  observed  pitch  rate  PFT  shows,  as  have  the 
runs  in  the  ether  t.wo  sea  conditions,  the  mysterious  double 
peak  in  pitch  rate  resonance. 


D.   SUHMABY  CB  BBEQUENCY  DOMAIN  ANALYSIS 


Comparison  cf  .measured  and  linear  system  theoretical 
frequency  response  shows  significant  correlation  under 
small-sea  conditions;  nulls  predicted  by  the  linear  system 
may  be  readily  observed  in  test  craft  data,  and  the 
frequency  demain  is  relatively  uncluttered  by  noiss.  The 
bulk  of  noise  in  the  system  is  not  "white"  noise,  but  tends 
to  cluster  about  sharply-defined  frequencies.  The 
techniques  developed  in  this  chapter  show  the  value  of  the 
frequency  dcaain  as  an  interesting  alternative  to  the  time 
domain  as  a  means  of  model  verification  and  validation. 
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Figure  34  -   SETUP  FOR  FREQUENCY  DOMAIN  DATA  REDUCTION 
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Figure    35    -       REPBESENTATI VE    WIDE-BAND    NOISE    SPECTRUM    OF    THE 

XR-3 
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Figure  36  - 


OBSERVED  PITCH  FREQUENCY  SPECTRUM, 
15  KTS,  AHEAD  SEAS 
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Figure  37  - 


OBSERVED  PITCH  RATE  FREQUENCY  SPECTRUM, 
15  KTS,  AHEAD  SEAS 
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LI-NEAfi    SYSTEM    PITCH    FREQUENCY    SPECTRUM, 
15    KTSr    AHEAD    SEAS 
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Figure    39    -      LINEAR    SYSTEM    PITCH    RATE    FREQUENCY    SPECTRUM, 

15    KTS,    AHEAD    SEAS 
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figure  40  -   OBSERVED  PITCH  FREQUENCY  SPECTRUM, 
15  KTS,  CROSS  SEAS 
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Figure  41  - 


OBSERVED  PITCH  RATE  FREQUENCY  SPECTRUM, 
15  KTS,  CROSS  SEAS 
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Figure    42   - 


LINEAR    SYSTEM    PITCH    FREQUENCY    SPECTRUM, 
15    KTS,    CROSS    SEAS 
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Figure    43    -       LINEAR    SYSTEM    PITCH    RATE    FREQUENCY     SPECTRUM, 

15    KTS,     CROSS    SEAS 
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Figure  44  -   OBSERVED  PITCH  FREQUENCY  SPECTRUM, 
15  KTS,  FOLLOWING  SEAS 
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VII.   SUMMARY 


A.   CONCLUSIONS  AND  OBSERVATIONS 


The  linear  tvo-degree-of-freedom  model  developed  in  this 
work  has  been  shown  to  be  a  reasonable  approximation  to  both 
the  Loads  and  Motions  program  and  to  actual  observed  craft 
behavior,  in  both  the  time  and  frequency  domain;  its 
relatively  uncluttered  simplicity  has  disclosed  some 
interesting  aspects  of  testcraft  behavior,  including  the 
nulls  and  nodes  that  characterize  the  craft  frequency 
response.  Furthermore,  this  simulation  is  extremely 
economical  of  CPU  time,  requiring  less  than  three  minutes  of 
CPU  time  to  generate  the  frequency  response  curves  fcr  up  to 
100  discrete  frequencies. 

The  non-linear  behavior  described  in  Chapter  V  is 
somewhat  supported  by  observations  of  testcraft  data  in  the 
frequency  domain  in  Chapter  VI.  This  interesting  behavior 
is  apparently  akin  to  that  of  a  relaxation  oscillator. 

Noise,  which  has  interfered  with  varying  degrees  of 
severity  in  almost  all  transient  analysis  done  with  the 
recorded  testcraft  data  at  NPS,  is  reduced  to  little 
significance  using  the  freguency  domain  method  applied  for 
testcraft  data  reduction  in  this  work.  The  major  source  of 
noise  in  the  XR-3  has  been  identified  as  primarily  ignition 
noise  from  the  eight  gasoline  engines  aboard  the  craft, 
transmitted  electrically  to  the  data  acquisition  system. 
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B.   RECOMMENTATIONS  FOR  FUTURE  STUDY 


The  linear  heave-pitch  model  developed  in  this  work  can 
readily  be  expanded  to  include  roll,  yaw  and  speed  effects. 
Completion  of  this  expansion  would  result  in  an  analytical 
six-degree-cf-f reedom  model,  and  lead  to  a  more  thorough 
understanding  of  CAB  craft  behavior.  An  example  of  the 
potential  of  this  approach  may  be  seen  in  ref  [9J,  in  which 
a  rudimentary  pitch  control  system,  based  upon  this  model 
and  utilizing  optimal  Ricatti  solutions,  promises  to  reduce 
sea-state  pitch  accelerations  to  less  than  ten  percent  of 
the  open  loop  response.  The  practical  problems  involved  in 
the  implementation  of  a  test-bed  ride  control  system  for  the 
X.R-3  testcraft  do  not  appear  to  be  severe,  aside  from  the 
physical  installation  of  a  suitable  set  of  stern-plane 
control  surfaces;  an  observer,  based  on  a  linear  model, 
could  be  utilized  as  an  optimal  observer  in  a  Kalman 
filtering  scheme  to  generate  some  of  the  system  states  which 
may  not  otherwise  be  suitable  for  feedback.  Since  this 
system  is  only  piecewise  linear,  the  microprocessor  appears 
as  the  most  feasible  means  of  implementing  this  control 
system,  storing  the  optimal  gains  for  several  operating 
speeds  in  a  table  look-up  or  interpolation  method.  Most  of 
the  system  development  could  be  inplemented  on  a 
slightly-modified  version  of  the  Loads  and  Motions  program. 
This  avenue  has  great  future  potential  for  expanding  the 
sea-state  envelope  and  improving  the  habitability  of  CAB 
craft. 

The  linear  freguency  domain  response  simulation  should 
prove  useful  in  the  study  of  the  non-linear  behavior 
described  in  Chapter  V.  It  is  felt  that  this  non-linearity 
can  be  modelled  in  the  complex  plane  utilizing  a  describing 
function.  Fortunately,  the  complex  gains  of  the  linear 
system   are  directly  available  from  the  linear  system  in  the 
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program  developed.   This  is  an  extremely  interesting  area  of 
future  study. 

In  the  area  of  parametric  analysis,  such  a  linear  model 
may  be  utilized  with  the  technique  of  root  locus  analysis. 
An  elementary  program  was  developed  for  this  purpose  and  is 
included  at  the  end  of  the  text  with  other  programs 
developed  for  this  investigation.  However,  time  and  space 
did  not  permit  its  inclusion  into  the  main  text. 
Furthermore,  such  parametric  analysis  was  inconsistent  with 
the  main  thrust  of  this  effort.  It  is  nevertheless  hoped 
that  the  program  will  be  of  benefit  in  future  studies. 

The  XR-3  testcraft  was  not  available  for  data  collection 
throughout  the  spring  of  1977.  The  frequency  domain  data 
reduction  techniques  developed  in  this  work  could  nor, 
therefore,  be  utilized  in  anything  approaching  an  in-depth 
study  of  craft  behavior.  One  area  in  which  these  frequency 
domain  techniques  may  prove  most  useful  would  be  the 
investigation  of  actual  testcraft  gains  in  heave  at  various 
frequencies,  to  ascertain  the  significance  of  wave-plenum 
chamber  interactions.  The  15  dE  discrepancy  between  the 
Loads  and  Motions  program,  which  assumes  a  plenum  wave,  and 
this  linear  model,  which  does  not,  indicates  this  to  be  a 
vital  factor  in  the  accuracy  of  any  model. 

A  serious  effort  at  ignition  filtering  is  certainly  in 
order.  Ignition  suppression  and  power-line  filtering  is  a 
relatively  straightforward  approach  that  can  greatly  enhance 
the  quality  of  data  collected  from  the  XR-3. 

An  additional  technique  for  data  reduction  to  the 
frequency  domain  may  become  available  in  the  very  near 
future.  The  MDS-8  microprocessor-based  data  acquisition 
system,  developed  by  Englehart  in  ref  £10J,  will  be  capable 
of  sampling  sequentially  and  storing  in  mass  memory  up  to  1U 
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analog  channels,  with  a  sampling  rate  up  to  7000  samples  per 
second.  Upon  completion  of  the  sampling  and  storage,  the 
MDS-8  can  perform  a  variety  of  data  manipulations  upon  the 
data  file,  including  an  FFT;  the  data  file,  both  in  the  raw 
form  and  after  manipulation,  is  available  to  the  NPS  IBM 
360/67  computer  installation  via  the  CP/CMS  time-sharing 
link.  The  possibilities  of  this  system  as  a  data  reduction 
system  are,  compared  with  techniques  previously  utilized  for 
data  reduction,  impressive.  However,  much  remains  to  be 
developed  in  this  area. 


148 


APPENDIX  A 


WIDE3AND  NOISE  DATA 


This  appendix  is  a  resume  of  the  more  significant 
wide-band  ncise  spectra  observed  in  the  FFT ' s  of  XR-3 
testcraft  data  collected  at  Lake  San  Antonio,  20  January 
1977. 

The  following  plots  are  FFT's  of  testcraft  data,  reduced 
as  described  in  Chapter  V,  for  a  maximum  upper  frequency  of 
200  Hz,  and  a  frequency  resolution  of  512  samples. 

STARBOARD  BNGINE  1/2  POWER, PORT  IDLE 

1.  Pitch 

2.  Pitch  Rate 

3.  Plenum  Pressure 

PORT  ENGINE  1/2  POWER, STARBOARD  IDLE 

4.  Pitch 

5.  Roll 

6.  Yaw  Rate 

BOTH  ENGINES  1/2  POWER 

7.  Pitch 

8.  Pitch  Rate 


149 


9.   Roll 

10.  Roll  Rate 

11.  Plenum  Pressure 

STABBOARD  PULL, PORT  1/2  POWER 

12.  Pitch 
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A     1    -         STARBOARD    ENGINE    1/2    POKES, PORT    IDLE 
PITCH    FREQUENCY    SPECTRUM    0-200    Hz 
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A    2    -         STABBOARD    ENGINE    1/2    POWER,PORT    IDLE 
EIICH    RATE    FREQUENCY    SPECTRUM    0-200    tfz 
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A    3    -         STARBOARD    ENGINE    1/2    POWER, PORT    IDLE 
£1ENUM    PRESSURE    FREQUENCY    SPECTRUM    0-200     Hz 
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A  5  -    PORT  ENGINE  1/2  POWER ,STA REOARD  IDLE 
ROLL  JREQUENCY  SPECTRUM  0-200  Hz 
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A  7  -    BOTH  ENGINES  1/2  POWER 
PITCH  FREQUENCY  SPECTRUM  0-200  Hz 
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A    9    -         BOTH    ENGINES     1/2    POWER 
ROLL    FREQUENCY    SPECTRUM    0-200    Hz 
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A10    -         BOTH    ENGINES     1/2    POWER 
ROLL    RATE    FREQUENCY    SPECTRUM    0-200    Hz 
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A11    -         BOTH    ENGINES    1/2    POWZF 
EIENUM    PRESSURE    FREQUENCY    SPECTRUM    0-200 


HZ 
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A12    -         STAREOARD    FULL, PORT    1/2 
PITCH    FREQUENCY    SPECTRUM    0-200    Hz 
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FREQUENCY  EOMAIN  LINEAR  SYSTEM  CSMP  III  COMPUTER  SIMULATION 


This  is  a  listing  of  the  CSMP-III  computer  simulation 
used,  with  minor  variations  in  details  such  as  frequency 
range,  tc  generate  all  frequency  response  curves  cf  the 
linear  system  used  throughout  this  work. 

The  CSMP  III  program  is  not  linked  automatically  to  the 
IMSL  Subroutine  Library  at  the  tf.B.  Church  Computer  Center, 
NPS,  Monterey.  To  acheive  this  linkage  and  access  the 
subroutine  LEQT2C,  the  following  Job  Control  Language  (JCL) 
must  be  inserted  at  the  beginning  of  the  program  deck: 

//SYSL  (3  DO  UN:t  =  3330,v:l=5E"s=0ISKC1  ,DI  SP  =  SHR  ,CS.N  =  SVS3  .  C  ^?3  .  IC-iCMOC 
//  CC  C I5P  =  ShS ,CSN  =  SYS  1  .FCRTL  IB 
//  DO  0ISP  =  SHR  ,DSN=SYS1  .MFSU  E 
//  CO  DISP=SHR,DSN*SYS2.iMSL.SP 
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PARAMETRIC  EIGENVALUE  CSMP  III  SOLUTION 


This  program  provides  a  root  locus  solution  tc  the 
two-degree-of- freedom  linear  model.  The  reserved  CSMP 
variable  TIME  is  redesignated  as  the  parameter  of  interest 
with  the  RENAME  statement,  with  the  range  of  variation 
specified  on  the  FINTIM  card.  One  additional  variable  may  be 
stepped,  using  the  PARAMETER  statement. 

The  parametric  eigenvalue  program  also  must  access  the 
IMSL  Subroutine  Lihrary;  see  the  Frequency  Domain  computer 
program  listing  for  the  appropriate  JCL  tc  acheive  this 
linkage  at  the  W.R.  Church  Computer  Facility. 


177 


LT 


CO 


r- 

II 

m 

z 

fN 

o 

o 

» 

o 

• 

• 

o 

ii 

<-> 

-a 

CM 

n 

II 

T 

32 

CC 

< 

o 

II 

UJ 

o  m  — ' 

•  •  >— 

<M  {^ 

II  II  h- 

i.J  l/>  O  *1 

X  ivj  • 

cc  •>  o  o 

•  •HO  • 

m  m  'jj   ♦  o 

CT  <JU  (_}(\J  II 

m  cm  x.  ii  — 

—  —i  ii  *a.  O 
'J><~                >r  2  u"*  — '  u 

s^-»            t  •■          'M  — <uj  "s: 

3U                "H  •                   t  ||  i/> 

•.■jj                ii  in               — I  ~>   -  ;— 

—  >              ««  m               II  o  "  «* 

CT>ij                a.  H                 -J  *U>  X                          v^ 

—  •-«  *.  ••  O  JC  — >r»-  —  O  I— 
•— <uj  au  u>»-*  ►  lmoj  •  ^- 
>—    »           ii     •  ►— O                 U>  ;\j.~ «  c/>                                                                          X            'M            U-J 

j_.            u^-u  *c  »            r-  •  '  "*                                                      u»        *v        a. 

3:— »!_)— «<tC?rn  *  II  «JJ  ^— »       "~             i— 

•  <j>          o.-«  1—    •               y>  ->0  -t                                                              t/3  -a.      —          ■*- 

—  >—          u_r\j  >ilt\                  •  »u-  <*                                                                       '_J      aj 

CT*— <                II  A,  II                 H Ol  ex.  x      *  '-J 

».U_            t/^«r  — )r*                  t»  »   •  LJ                                                                                  -    t         1              <^ 

O"-*          •—  □»(_;►  _<  f>    •  2:                                                                  !fi>*      O          < 

—.  <           X*  O                 »  "^O  —.—.-<                                                                  x—      ~ 

>-«    »            ,—  vj-  ,—     •                      •  •  f\J  O                                                       UQ                                                                                          CI.  l^J        <-£             C/1 

rO-«                 •  U.*3               "3  I    0  1    •                                           "-I   tX                                                                       *       x.             ij 

»CT             ««C  — «  <T'>J                   OC  ■**0»  1U\                                                       'A  £»—                                                                                     v^  w         O                j 

-»—          -.  11  'x.  11                11  <n      11   11  iu  11                                       pjou                                                           □*■       •           -£. 

0s— 1              <  UJ              oj  a      a>  hoi                                        ujuj>-<                                                                   ~t      <T          U 

•  ^           OS  ►                 £  UJ      f_)>  UJ                     2T                          C5d^:                                                                       U>»       I             U, 
^  >           u»Z  u.:0>0|—      X  —«  5!    ••  "3                -»       *  #•                                                                           0_jr\l 

—  O           -">  ^l  _!•«         »LU        •••■  ~xro  ■-*                >-«      — 4<N_i                                                                      — )<       *■  <I      'LJ 

,jj              •  <;— 4      u       •  <i     00  5  •            <r             *jcq<                                                      :>~      £^     O 

-r                   '/)0  •-•  II         2-        O  T.        Om  CJ.nO                   INI                    UlU  — >-i                                                                                           LJ       — '  3f         ^>        31 

q  x        i—   •  cc  x     *-;     >o  <i     u~\  •  n             1—              »QZZUJ                                                      <  —     k  «^     ^     < 

>-<ui         ZM  «jv—     uj      11  a.      0  H  zu.             _j             rvj^-rf  j-c                                         oo             w      ai-n     iu     03 

i/o-J          ^m  ^»Q     cj     — i  ii  _j  »— « __j             <i              wtti_Krt                                             •   •         i— +  accjcr      .£     ai 

^ira.        ;— ii  i  »—            ct  i—     3<  u_o            -<            xsssi/i                   j               c— •        jo<i/^.     a     a. 

•LUi            LOO  ^TiJCJU.  H-  •0_l_l                                i.       C^O^Oll    II             iJJi-i  CO   I    •♦■                  ■*■ 

~32:o         -Z  7)          _J  <      irar  _jct             »-«              o   »^j:uj                      -—       »  ••  •-»—         _j30j^oo      >aj     oj 

•■>—  o           OH-  "Z>—       <      !—  ^      XIXJ  3iU                2T                nO— *  I     I  O                          I—       -i^iOT--            j  —  O.  I      •      i—        i 

>— a     _j    or:  ^     -£     ^  t_>    k_j—  jc^           ^_j        r>\  » >sn/^<i     -j        o-*      h  u  ii    -  -        >-r>*  » ~--^     o      i 

<         «i  <i     uj     <i  lliUj  i— uj                <»—      iitj-Siix      #         ►—^i_^_-).^H^t_            oo  ii  ii      a.      ii 

O          >-•          h-  ►—      Q      I—  SX  2T2.                    i-tQC      0<I  II   II  -u      i/io      2!-C  £           -3>-'^'                 •cx.ri.-~i      X      <S> 

uj           j—           oi  r/>fc_ioo  «*<'C:*:T  H-O      ujkOO>       ZS      <   •CtT'1^    *•— •— '                fN«K<      U      '-" 

^             »-*             Z  *Z.        I/}       4.  CCUC  l-JLT                         >-<00        C_10— 'CnJ<1         II   Jt       Z— "VJ             •— >i)GD                     I    HJXa.        O        nX. 

-<        2:        c  vZ           a  <i<  <t               zj     <  utoui        coii     >  ii  -oo—                  ii  ua:cr!           a 
u_>v»>».-«#*u-»#*o**  *t_)*»  *  a.  a.*  »q.»»-»  #  «-«^-».cic^-S3:»  *  <i#a>zuQ<ir\u<*  ^>i-  aca.Q.a.*  #  *x 


178 


vC 

XI 

cc 

< 

UJ 

o 

cc 
a 


to 


X 


XI 


— » 


_l 

cr 

I_l 

* 

(_i 

«£ 

U: 

•=4. 

Jk 

'X 

> 

UI 

I— 

a 

X 

X 

->:■ 

o. 

Q 

;c 

a. 

X 

_i 

C_S 

rr 

►—* 

* 

— 

a. 

VJ 

_> 

a 

• 

-j 

<NJ 

.j 

Ui 

«_ 

n 

UJ 

a 
.u 


a. 

X 

•to 

LO 

UJ 

ox     >— 

M    * 

<35 

to--. 

UJ  II 


o 

m 

•      LU 

H       Z 


o    z 
— ■    n 

II 


UJ 

to 

z. 
cc 

UJ 

V- 

(-0 

XI 

a: 


UJ 

(/; 

3 

a 

a: 


UJ 


Q 


•u 

X 


^1 


O 
(J 

t 


1       * 
UJ       — 


x> 


r\j 


r\j 


i 


u 
I 


u  u 

x  <yi 

•to  x 

•      — *  •        ♦  ^j  + 

(Mr-t—JtN        "\I<"M_;     •  f\J<\| 
>*_»<>*      >*_J<ttVJ_J-J 

_J<UJ_J  ..  J^UJN^^. 
>*UJ(^— »_J  —  UJ  O—J.iJUJ 

i  wm*  «*  *  (/i<j.  — »ooi/i 

rr:  _i_J*cr_Joo_j-»-Q.'>j+* 
"3«.— X<<t      —  X_J— \-*-* 

i/)OQaj(/)'KOQtU'jjt(< 
UJ3J  »-<Q.IUJ  •  *-»«/>  l/JUJUJ 
KON3  I  I  I-AJ3:  1  X'^-O 
_1  UJ  ||  ||  n  ||  u~)  II  II  II  II  J.LV. 
Ql  —(—.^-i—j  (\jrNJ(\J(\J  II  II 
"Z.  CC  ~J  _J  — I  _J  CC  _J  _|  _l  _l  _J  _J 

ou<«<jn  <<-*<<< 

(_5  U_  UJ  UJ  UJ  UJ  U-  UJ  UJ  UJ  UJ  UJ  UJ 


X 

re 

i_ i 

□C 

<-J 


c_j 

L'J 


uc 
~) 

> 
x- 

H- 
UJ 

2: 


t 

JE 

XI 
a 

u> 


•to 


LJ 

*CC 
'^.CJ 

.—  * 


3< 

o 

_J 
* 
o 

f\J 


<i 
•X" 

CJ 
jj 

LJ 
•* 

—4 

ix 

l_J 


<\J 


x  c 

■■—  aCDJS. 
*  —    •   I 

ex  II  3: 

(JI-pH  II 

X*  •■^— i 

I   0_5  .O 

CNJ     .~^ 

■^f\JO  — 
-J  —     •!/> 

i  — *  •  • 

-K-JO 

«  a  •  i 

t  •  i-^G  *-*•—* 
jEii-Ji/OcO 
C«I  II3S 

OtOXX 


CC 

X 

(— 


'\J 
1/1 


< 

r\j'^> 


rsj 
i: 

<i 
aj 

LJ 

_J 

<\J 


a 
< 

:jj 

* 

re 


-»■        * 

•UIVJ        * 

Xl-x      to 

.— CJO-* 
■i  —    •    I 

—  ^  O  UJ 

3<  ii  3: 

_}V_(\J  II 
X  *  .•")  f\J 
•to  O-S^' 

V.INO  — 
_J—    »Ui 

—  ^O  JX 
■to  fNJ     •     • 

J<JO 

II  OJ     •     • 

CNJQrvj  t\i 
X"_l^^) 

<r  H  jus 

CCfNJ  —  — 
Ovoxu. 


2. 

a 
^: 

<i 


o 
cc 

C 


X) 
CC 


CVJtM 


<  -I 

I'JUJ 

xx 


7*2: 

to-!— 
*    * 

eJ<^» 

XX 

I  + 

•  • 
>-NJ<M 
^\ 
-J -J 


ou 
x-«- 

I  +-• 

•  »>-r\J 
•NJtM  -  > 
W_J<_' 

_l_l  23 

*  *  *  X 

-HtM  —  -H 

33oO 

*  -i-  xo 
OC  •  X 

*  •»    •  + 

f^(_JCNJ     • 

i-X^fM 
CC  CC  _l>v 
«  »  — '—I 

oo*  •• 

•  »ir\-» 
cnj.-m  »m 

lit* 
H  ii  ii  it 

— lt\J^<Nj 

'j-;^u 

3J3J 

£Tjd:cG 


a 
o 

UJ 

o 

JC 


< 


u 


(M 


* 

J> 


X 

a. 

■*■  s- 
—  <J. 
o_i 

•  3. 
(MX 
\« 

>o 
o  • 

I 

l|  II 
2T2 

•rf^ 
X_J 
XX 


o 


(NJ 

•  T 
(NJ 


y 


CC  o 

I-  I 

X  M 

o  ^ 


Lj  ii  y^tot/JwO      i/it/iio^ooooT  Otoxx  Ovoxu.  i^orcco  xx  __ 

*0^^***#x<xx*x<fxaxx#*»**_i^i--ii-i**#x3:>-i»--i#»*xxxx#-»#xx**-»»x 


179 


UJ 


2 

O 

O 

CO 

UJ 

_; 

3L 

J£ 

O 

JJ 

a 

UJ 

t^ 

000 

3T 

m 

X 

X 

-12 

n 

a 

L- 

1— < 

z 

►- 

•  .—  <c 

< 

<r 

r- 
< 

2<a  u 

i/> 

i— 

•-to-S 

UJ 

UJ 

i 

j_2      ~2 

o 

X 

3 

Xuj   «■ 

x 

_ 

h- 

^: 

X< 

,— LU 

UJ 

O  —  •— 

X> 

C! 

—1 

^:     uj 

o 

•■* 

2 

'X 

— 4UJX 

Li. 

-JUt. 

=? 

•—»—.— 

21 

^1— I 

•« 

•  • 

UJ 

•-^_;  » 

-5-f 

3 

>- 

< 

h- 

I 

(XU.L— 

♦  if 

CJLJ 

u. 

t- 

X 

►— 

"JiO 

2  -■ 

— </— 

<£ 

•u 

•--4. 

a.  -33 

1/1 

-i/ 

X 

^- 

2T 

3UN 

XC. 

UJflC 

6 

L- 

3 

—4 

-    •> 

X* 

XI— 

X^ivl 

+ 

UJ 

CJ 

3 

u. 

QC 

X  -J   » 

£.*£. 

U_X 

*— 

>— 

a 

^—4 

x— t-n 

<*.•«* 

u>- 

pM 

<1 

1—  UJ 

-J-* 

i— 

— 

O') 

r--U 

U.«— _J 

Q.U. 

'•ii_) 

<_> 

a 

•jj 

l»o 

X 

'JUJJ 

XU- 

OJ^ 

ai 

UJ 

3 

(_J^v 

'J2 

■— ■ 

*£*X 

»  •» 

->•*<. 

u. 

a. 

—I 

>.«* 

■~ ■ 

*;-_>  — 

•  • 

ivrvj 

>— 1 

on 

CJ 

<4 

--» 

V) 

<r 

'JOC 

CO 

>■*• 

►— »— 

LU 

UJ 

> 

_jr\J— 4 

o> 

a: 

•-.     <», 

!— 

— >»_) 

<r.x 

X 

T 

GWi" 

-.r 

— ■  T" 

^-?-^> 

;?: 

33 

XC 

15 

|—        -"NJtiJIU 

2 

-»   i 

-3 

'u 

rrjaj 

-•at 

X 

X 

— i 

'JJ 

uj     Q(f  at  ai 

<e> 

JT*-iUJ 

2. 

ia 

XO 

1— 

L_ 

u 

*H 

IhU*  * 

j 

>->< 

»—  Cf.  >— 

J 

*  + 

JJ 

-4 

-4 

-3 

o 

!  -    U  ">  X  _J 

«t 

UO 

X  -0    i 

^. 

—  — i 

— >  i.i- 

S. 

s. 

_j 

• 

<t. 

uxx>-u 

—4 

j_   » 

*IW»- 

>>- 

U 

—i 

OLJ 

II 

J 

V_>—  '-L_jOJl~ 

- 

— -  w^ 

UI_Jc/) 

C_J 

CU 

—J 

IS) 

o 

oo 

^_l 

|— 

u. 

" ■  C3  UJ  -^  t  Uj 

»— 1 

l1-« 

t— 4 

2 

33 

<t/) 

UJ 

o 

X 

—a 

UJ 

O^--  ^TX 

2 

T»* 

lU^JJ 

< 

crrc 

-»LU~ 

■> 

• 

> 

Ov. 

X 

r> 

r-l «—  Qt"  a  +  H- 

1— 

<* 

X3X 

xa. 

■—2D  rr. 

t— t 

II 

*— 4 

r^— ► 

r~ 

-; 

00OC  +  -JO 

v—  <-« 

l-LUJ- 

oo 

♦  * 

x— jiu 

1— 

<_> 

V— 

-y)0 

□ 

; — 

X 

'JJ  fM   1   X  Q  * 

UJ 

#  <l 

UJ 

_i_) 

<'i"Nj 

M. 

-J 

-t 

LUCNJ 

— 

LU 

h-    • 

(£i/>XH-— 'm 

X 

J5X 

-  JC 

u 

«-!.  <1 

X  J> 

> 

C_J 

> 

Ulf/1«*.tJX 

'JUl/I 

1    IU>— (— J.^'O 

;— 

(JUJ 

LUM> 

X 

UJ'ij 

1 — 1 

*: 

►— 4 

— » 

—4 

1    JJH- 

• 

1— 

s:— j 

— i^O'NiQ    • 

o 

X<2- 

*—  X 

a 

00 1/> 

tUJfct— 

» 

<X 

o 

ex: 

— ITCOJOQ 

< 

t/?  I  <\i<r«  —> 

_j 

u 

+  ~" 

X-<L_ 

u_ 

ia 

XUJ 

— . 

lU 

a  • 

LU 

U)   |    X 

• 

+ 

^_j 

U-  CN  (i  r-3  f\l  "• 

CJ 

rmm 

u* 

JHU 

♦  + 

t-2  o 

^T 

o 

-JO 

c 

'JJOOI— 

— > 

vl 

-  a 

•X  f)  O  -It  «/1  < 

H-  X 

3 

_jro 

a.  vOUj 

_ i 

</)<•> 

_i 

ft 

•*     •      — 4(N 

x.oo# 

LU 

;_^(\j 

j:— • 

—  'iJ-rt-  — 4'Xli— 

L03 

X 

— -> 

X       X 

'X 

UJ'.iJ 

UJ.JJ«T 

rn 

_l 

-^<-»QoOOO 

_J 

— .UJ— < 

• 

UJOOoO 

o-i 

II  iX—HO  YTUJ0D2: 

r 

Vc  * 

U'  lu) 

3 

aia: 

1-13 

* 

<t 

<~>X_JliJUJ 

< 

ii  j:o 

u_ 

Xujuj 

•—  'JJ 

Xw^X | X 

•  <— 4 

ai< 

U'-'  JJ 

— i 

a.  a. 

3t-oj 

rvj 

►-• 

•   •car.  a:      *-> 

•.-H 

_J-*      1LU 

i—  cm:     — > 

ooor 

t—     II    —"(X^L,. 

r  *|'« 

o 

<u 

XX 

— 

XCL 

U.        — « 

— 

L- 

ii  a  +  ii  ii 

X 

Q  11    II 

X 

ii  H  ii 

(X 

OXa:  i    ii  il 

4-  r" 

1  xo 

L/)|— 

OO 

II    II 

ro'/i 

a 

aC 

Qjaoo    u 

a: 

— <-lh- 

1— 

<oo    c 

i—LU 

— if-  a  wv_<  o  r- 

-"XCC 

j— •       X 

u 

-4  CM 

o2uj 

»— 

<r 

JQJ- «r\J       H- 

■<* 

OClJ'JUOH-—"^       ►— 

<TX 

arc  ii  H  Luf— 

JO 

>  II  X 

X  f> 

JC 

KT^i 

u<rx 

Ql 

(_3—  II  OOOCNJ 

a. 

QOsIX 

■* 

ujuooo-n 

tu^- 

LJ<NJ»— ox-JJ 

II 

II  O  ii 

t—  ^1>— 

XlU 

r^ 

LL-OUJUJ   II  lJ 

XV- 

XXujuj  II  "3 

•X  UJ  _l>—  X 

30 

tco-a 

U-_C 

Lk  ' — IIU  '_,  LL.  LJUJUJ    II    L_ I  i'-'i.XUUUJ     II      J  X  JJ_J(—  X  1  <— '  iUJ^J  L^..» 


180 


t— 

< 

a 

x 

c. 

5:' 

CD 

<t 

2: 

0 

| 

Q 

1 

2 

# 

<t 

** 

^^ 

m 

< 

< 

h- 

LL 

u 

_; 

j. 

< 

;— 

m 

# 

1— 

e> 

r- -i 

'j 

CJ 

X 

^M 

+ 

INI 

0 

0t 

J 

a£ 

— 

<X 

> 

33 

_Q 

32 

«r 

a. 

t 

a 


u 
at 


< 


O 

OJ 

X 


I* 
< 

CO 
CO 
Q. 


a 

♦  * 


an 
a 

LL 

#  *  * 


OJ 


>  + 

*-«.'.TJ  "X 

3:— o. 

<»  + 

coc;.x 
coo  "5 

axes 
—  *  er 
^  id  a. 

OJ<r  — 

**■*  * 
Ly  II  II 

II  IX) 

•nih-x     as 
airoaa      O 

an     u_ 

*  QCO* *  *  * 


5: 

OJ 


t— UJ 

fMJ_ 
»VI 


j-oj  nc 

h-3C  UJ 

iniCOCO  ►— 

X:LLU. 

O-OLj  > 

a*  *  a 

♦  cunj  3 

xx«r  xi 

oil  _l 

I  11   11 

II  ICO 

•VIS-  3L 

xxx 
000*  *  * 


<Nj 

— >  «• 
o  — 
00 

x<_; 

+  x 

•  - 

fsl   • 


i/;r\) 

_s</"» 

*  _S 

tyM 
X  * 

•  o 

>»  >< 
_j  1 

-— »   • 

*  CM 

-*v. 
C/)_J 

3S<-» 
«»# 

a-  -12: 

*  -idJ 

x* 

CCCJ_J 

»  *  <t 

OQIJJ 

•  I") 
CM  (XX. 

I  » 

II  X 

•VI  — 

COX 

XX 

*  00*  # 


t3 
O 

I 


< 

13 

f^CC 

ri* 

CD -J 
O-X 

*  I— 

— <LJ 

LJ(_J 
X* 

UI 

*■  f— 
-1.0 
>-■■— 
!-!.£ 
•-«  I 
-£  — 
I  — 

■ni-3 


<1 

OJ 


LU 
'./) 

Q 
Z 
<i 

a 

CO 


— »LL  Q 

rj-*  — ' 

X— t_J 
•    -»< 

rvjujo: 

JX* 

—  ♦  X 
■*  Ik— 

o>— c 

•  Qh 
•OCM— «JE 
<i  (Nl  •  «  "5  I 
U>        II  <nj  II    I    II 

00  — «    II    — •    II    r-4 

-5     J-*I-»I 

13      <i>j'_inji~ 
cu      lu  q  t_j  v/>  '.o 

wxxn 

»  *  *  XOUQQ- 


•■NJ 


I 

a 


o 


UJ<-J 

^r  11  osl 

QC  CNJ  II 

xi_jr\i 

VJLULJ 

*  xc 


_!"-■ 

<*.t  "> 

— .  UJ# 

O-^  X 
xr\jr— 
+  |vit-J 

•  OX 

(\J>C  '_) 

^O* 
— '■»  X 
—  -LH- 
■*  S-C^ 
OCJ>- 

•  ►-  JC 

AJ3  I 
II   I   II 

fM  II  (\l 

INI 
I «J<— 

XII 


O 


a. 

a 

X  _J 

:—  --i 

CM  t0  UJ 

NX'-" 
C-"_3X 
X  +  * 

a-<x 

+  X'— 

ivlt.0>-« 
i/lU 
xj  l 
O  II    II 

II  Xro 

vjh-5: 

III 

(-103  * 


to 


a, 

»— ' 

< 

_l 
a. 


* 


CM 

'  1 


— 1 


CjC 
fM 

> 

o 

I 
II 

I 


r 


o 

X 

» 


.1 

J 

<_> 

a: 


^1 


*»**I3*** 


c       : 
rsjrv 

^r  V. 

*  W. 

>*■   ■ 
**■   -» 

ia 

at! 

*  zc 

LL* 
—ILL. 
O— I 

•  o 

II  I 

I II 
I-  < 

a.  ti- 
ll 
3C 


XI 


3 


(_:,{.-* 


181 


h-S 

O.XX 

iQC 

5:c*  ■* 

nc*  -»— ' 

OXX 

_u  |  UO 
a:ooo 

O.OXX 

ax— «• 
— ■  n  ii 

II  TCD 

rsjt-  — 

UtCLU. 

a.  a.  a. 


rn 

— >* 

(\J-S- 

*  «-• 
*o 

—  o 
ox 
o  + 

X     • 

+  -\J 

fM_J 

wi^ 

*  Ol 

C\J-X 

V}  •* 

-SIM 

I    * 

*  — 

•»■  o 
— <J 
ox 

O   I 

x   • 

i  »>* 

rvl I 


*  ;/) 


LU-2  * 
H-wO 

O*  X 

DI   I 

O-Tf  It 

II  J- 

Nt- 

COCO 

ax 


< 

lo_i 

a  o. 

xi 

o 

CO 


X*r 

c-« 
*  _l 
— <  -I 

<lO 

LLiX 

x  — 

+  w 
— ■'_> 
OX 
O    I 
X     • 

I  "N 

-Nl— I 
^  _l 

_l   '-I 

.ux 
o-  I 
X  • 
I  CM 
•>s 
r\j_i 


lu         w. 


•»  X 
uOX 

xo 

O  I 

I  II 

II  -* 

INI>— 

a.  a. 


LU 


uj 
'— * 
x 

i 
o 


fNitN         — < 


N   t 

LUX 

x   • 

+  <N| 


j_X^ 
xacc 
v.^X 

oca:  £ 

23  13  X 
COCO< 
XQ.CL) 
*  #  -U 
OUU. 
•     •■«• 


X 
x 
u 


>• 


X 

u. 
to 


1 

iu 
/) 
X 

I 
r\jo 
X    • 
•  •Sv       r-<N 

— »<NJX  •*  -» 

*  X  t-J  —  oj 

cni— +  ii: 

N/)"i(-'0 
'/I1NWX 

O  II  <-lc_1  J 
II  :\IQ  II  n 
(NX  II  XTO 


•u 
a 


o 


Q. 

3 


n 

x 


a. 

a. 


i_ 
5T 
LU 
^r 

i. 

Z 
•— t 

X 

_l 

cm 

\- 

o 


~fc-  -X.  lx,  -L    —  -u  -X.  <-x-  :_  _  i_w  t-i_  Ihl,  -t.  UU 

*  uqo4 *ao***#co**#  qoqoo*  *  *  a#  *  *  * 


f— r-(>J  .—  — » 

iX^w  —I  X 

l»oi-  a.  x 

./)/■)  ^  x  >•    a 

X                      \^'-*  o  >      + 

u_                         jj/)  *  '-•X 

u..*:£            <i<i>v  2.     x  "s     — 

X  i  X             f  ->.ji  \     ■—  — •      i     > 

<»■  c.lj            oo-»x  »i     x     2:  •ji     o.     >-< 

—  *  *             ++«£cl  i^i     o     \  a.      +     >■ 

o— —            j-foa  x     •*•     —  o     x     — 

0<J)0                 C-'-"*'    ■  C       X       3D  ■*■       h-       CZ 

Xu<-J                IKI-  ■»•       !—      21  <*sl.>-> 

■♦■■^X                  OCoiOl— "^  "^       CO       3_  >        CO       "^> 

0  +  +                 ;*Ix    •  CO      XX  OO.Q. 

•  uo           ^^rix^n  x     '_>    o  u.-»-^>-o 

— t    •    •                — ,-4..2^  C       *       *  0-  +  >  + 

I  — <  — <                 "**.<      lu  +I11J5  *•— X"-"0 

O    I     (                     ^-Q— »f\J  i-M>s,j_'s.3;  'X'V^'vi.        XJ2 

•  oo  ^5:*^—  asaayi  o  -£o.o.a*^i—  >. 
<m  •  •  j_     cocc^oe—  ooiniiiiooaiiaLLaiitta 

'v-rvjfM  r-      o.a.03— »<J.  •    »OX(OXO    •    lOQ-Qaoaaa 

—  _j_J  CO    I  I  C^'i  II  II  II  II  II  II  II  It  II  II  II  II  II  II  II  "  i> 

II  w—  2!  II  M  I  <— C  -.-*—.-.—  —— .—  —  -^  —  — .  —  -.-— — 
XH  it  X33HXIX  fM-n^jr^.^^o^^orvjrou^OL^  ^M^-J^ 
>-j:cl  ^     f—j-tvii—  »•»•••»«»•••••••• 

lululu  a      coc0co>-i— <  — « r\i  m  rn  (■*>  co  <"<"(  vf  u>ohJ^J  v0^or^>-r- 

lai  x     sj-'lYO  —ww— v—— ^.— www^— — — 


182 


UJ 


X 

-S 


tr>  -J 

►  u 

o  c 

UJ       — >  •• 

>     -  o 

UJ      ..  • 

>^       UJ  — *  >■ 

UJ       UC.  "•  >-  :_} 

■       M.  -Jh-  <£ 

_i  uj>-«  uj 

>-     oo  a>  x> 

jj       oj  r— »—  u 

—                      2J       3  JH  UJ 

•>     -J  i_j-«  <3-                        >■ 

•  •                  — <     <                                                ►</)  j_  o 

x                 u_     >  <->^  2: 

— .                _i     «r                       — •     ^j          »lu  o  _i                   lu 

a:     — 1          <i      lu     —                  ^     ^  cut/-,  «£          <t                      3 

t—     o*           •     'J5     -h                   S5ii  -«  <                       J                      :—          X 

<!»<_)'— «     -•                   ■>•»>*  uj  >•  a.          _3                       ,,J          — 1          **■          *— ' 

2.        — t                •       UJ       — I                         — »       ^— ■  >_H-  2.             1—                                QC              <3.              <l              t— 

II                 0s                      L_         — .               •*       *»—  ►_•-«  <3                1*                                     U_                1—                                     V- 1 

<     -3         *     jj    -j—-         n    in«t  1—  o  a        zr                            h-        Q        a. 

ijj        — 1 .— »          ••       t—       ««v— -         — — . 1  — - » I          «•_!  □             Q  2T             «3             LU             UJ             LU 

X      ~)in      :-.     -       —>_<••     r\j <r in<: --i.  ouj  c          O          *"*          cl          T.          :—          >— 

u_       »  •     cC       •     .-•  «i  in      —  >-<;>>        .  >  -<         w  q.          3          3          -j.          <* 

"-rn       •     x     — ».>  •      _j"?_ioo  '-J  x           1:          x          1—          ^-          i—          ■— 

•     — .-*     o»     — )     -join      -i^<h-v-<  <^ui  uj         uj         <         <i         x         uo         y> 

>«.     <u_       ••       »     u..-.-*     >uj>'jjuj      ||  a£  a  a.         a          ."2         'jj          1 

Xco         *v     — v*      <t     ^      _jujul     O"—  O— —  2.>— •  f\j                      r*l                  >-<.      >- 

ivj»—5:     — »«—c^-—  tjjw- -\cr<-.c\juj>-<_j -•_)_)  -«►-.<  >—        t— <j     x        x<j     t—        o»—     c 

oajm-o     o       «nj     Oil^-  ••  11  in— zjuj-T'j-K-3.  ^-.ZK-ujnc— «t_ui£i—  ujofNJujij>— jju.     uj  <r  uj  uj  <t 

•  2.2.2:     >r!— -»sr»— z^jstk--^— >rt-2:— uj'— -.'juj  ^--iijo'jzoaiijcsi-zu1-  u.'0<qoiuio^ 

nQQQ         »<C  II     »<MO    •»<  II  — <    *< -<cOa:c"),jCa.  -"ON^i-  -S^Xrvi  ji^Jla:    -«^UiLr._lDLf—  *—  at  — 

II  II   11    11       sJI'-'-uii-w-OJ.HH—HJSi-ti  1  cu  1    1  u- 2.      u_i/>      ido      tua      uju-      iucj      aji/)      uj«/j 
— — — — w-— cr:     —otr^f  tij—a.     j«»a^<iii«iii  11  ^5:     >-2"     i— 27     ►-2:     h-2:     1—2:     t—X 

o>  nj  in  cr  cc  iu  -3  — «uj  do     u-iao-a  ■ uq  du-^u  in<"<">  ai— ixj  _j_j     _iZ)     _j-J     -jXi     _j_) 


183 


LIST  OF  REFERENCES 


1.  Booth,  F.B.j  The  Frequency.  Response  and  Operating 
Characteristics  of  the  XR-£  Loads  and  MQtions  Program, 
M.S.  Thesis,  U'.  S .  Naval  Postgraduate  School,  Monterey, 
Califorria,   1976. 

2.  Leo, E.G.,  and  Boncal,  R.,  XRzl  Surface  Effects  Ship 
Test  Crafty  k  Mathematical  Model  and  Simulation 
Program  with  Verification,  M.S.  Thesis,  U.S.  Naval 
Postgraduate  School, Monterey,  California,  December 
1973. 

3.  Menzel,  Rein  hard  Fritz,  Study  of  the  Roll  and  Pitch 
Transients  in  Calm  Water  Using  the  Simulated 
Performance  of  the  XR^3  Surface  Effects  Shi£  Loads  and 
Motion  Program,  M.S.  Thesis,  U.S. Naval  Postgraduate 
School, Monterey,  California,   December  1975. 

4.  Mitchell,  H.  R.  ,   Investigations  of  Methods  to  Optimize 
.Captured   Air    Bubble   Surface   Ef_fec.t   Jhio   Digital 

Simulation  for  Irregular  Sea  Conditions,  M.S.  Thesis, 
U.S. Naval  Postgraduate  School, Monterey,  California, 
September  1974. 

5.  Gerba,  A.,  Jr..,  and  Thaler,  G.J.,  Pressure  Ratio 
Effects  on  the  Heave  Motion  Characteristics  and 
Pressure  Dynamics  of  the  XH-:  '.  Loads  and  Motion  Prograx 
for  Stec_  Weight  Transients,  Progress  Report  to  SZSPO, 
U.S.  Naval  Postgraduate  School,  Monterey,  California, 
January  1977. 

6.  Speckhart,  F.H..  and  Green,  A. I.,  A  Guide  to  Using 
CSMP-the   Continuous   Systeja   Modeling   Program,     pp. 


184 


220-235,   Prentice-Hall,   1976. 

7.  Gerald,  C.F.,  Applied  Numerical  Analysis,  p.  16, 
California  State  Polytechnic  College,  San  Lais  Ctispo, 
Califoriia,   November  1973. 

8.  (No  Author  Listed),  Installation  and  Operating 
Instructions  for  the  TDC  SuperPanel-50  Spectrum 
Analyzer,  Instruction  Manual,  Time  Data  Corporation, 
1976. 

9.  Mclntyre,  L.F.,  Optimal  Pitch  Control  of  the  XR-3 
£ap_tured  Air  Bubble  iCAB}_  Test  craft,  Private 
Correspondence  with  Dr.  D. H.Collins, PhD. ,  U.S.  Naval 
Postgraduate  School,  Monterey,  California,   June  1977. 

10.  Sngiehart,  CD.,  Da^ta  Acquisition  System  for  Unsteady 
Aerodynamic  Investigations,  M-.  S.  Thesis,  U.S.  Naval 
Postgraduate  School,  Monterey,  California,   June  1977. 


185 


INITIAL  DISTRIBUTION  LIST 


Nc.  Copies 

1.  Defense  Documentation  Center  2 
Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  0142  2 
Naval  Postgraduate  School 

Monterey, California  93940 

3.  Department  Chairman,  Code  67  1 
Department  of  Aeronautics 

Naval  Postgraduate  Schaol 
Monterey,  California  93940 

4.  Department  Chairman,  Code  62  1 
Department  of  Electrical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93940 

5.  Professor  Alex  Gerba,  Jr.,  Code  52Gz  5 
Department  of  Electrical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93940 

6.  Professor  George  J-  Thaler,  Code  52Tr  1 
Department  of  Electrical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93940 


186 


Asst  Professor  D.M.  Layton,  Code  67Ln 
Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California  93940 

Mr.  A.  fi.  Anderson 

PMS  304-31A-1 

Surface  Effects  Ships  Project  Office 

PO  Eox  34401 

Bethesda,    Maryland  20034 

LT    Lewis   F.    Mclntyre,    rJSN 

21    Burobough  Place 

West   Asheville,    North   Carolina    28806 


187 


20  FEB79 

3     JUL  7  9 


25210 


Thesis 
Ml  8955 
c.l 


171209 

Mclntyre 

Investigation  of  a 
linear,  two-degree- 
of -freedom  simulation 
for  the  XR-3  Captured 
Air  Bubble  (CAB)  crafi 
in  the  frequency  do- 
main. 


I 


rEB79           252  10  )9 

i 

■ :: 

Thesis 

171209 

Ml  8955 

Mclntyre 

c.l 

Investigation  of  a 

linear,  two-degree- 

of -freedom  simulation 

for  the  XR-3  Captured 

Air  Bubble  (CAB)  craft 

in  the  frequency  do- 

main. 

thesM  18955 

Investigation  of  a  linear,  two-degree-of 

I  II II II  NIMH   II 


3  2768  001  88192  3 

DUDLEY  KNOX  LIBRARY 


